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ABSTRACT 

When a jet is perturbed by a periodic excitation of suitable frequency, a large-scale 
coherent structure develops and grows in amplitude as it propagates downstream. The 
structure eventually rolls up into vortices at some downstream location . JVe^pproximate 
the “wavy flow” associated with the roll-up of a coherent structure by a parallel mean flow 
and a small, spatially periodic, axisymmetric wave whose phase velocity and mode shape 
are given by classical (primary) stability theory. The periodic wave acts as a parametric 
excitation in the differential equations governing the secondary instability of a subhar- 
monic disturbance. 

The (reson an t) conditions for which the periodic flow can strongly destabilize a sub- 
harmonic disturbance are derived. When the resonant conditions are met, the periodic 
wave plays a catalytic role to enhance the growth rate of the subharmonic. The stabil- 
ity characteristics of the subharmonic disturbance, as a function of jet Mach number, jet 
heating, mode number and the amplitude of the periodic wave, are studied via a secondary 
instability analysis using two independent but complementary methods : (i) method of 
multiple scales and (ii) normal mode analysis. We found that the growth rates of the 
subharmonic waves with azimuthal numbers ft _ = 0 and /3 = are enhanced strongly, 
but comparably, when the amplitude of the periodic wave is increased. Furthermore, 
compressibility at subsonic Mach numbers has a moderate stabilizing influence on the 
subharmonic instability modes. Heating suppresses moderately the subharmonic growth 
rate of an axisymmetric mode, and it reduces more significantly the corresponding growth 
rate for the first spinning mode (i.e, fi = 1). Our calculations also indicate that while the 



15 


presence of a finite-amplitude periodic wave enhances the growth rates of subharmonic 
instability modes, it minimally distorts the mode shapes of the subharmonic waves. 
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Chapter 1 

INTRODUCTION 

1.1 Motivation fa Technical Importance 

It is now generally accepted that large-scale coherent structures, in addition to 
fine-scale chaotic motions, exist within all turbulent free shear flows. Free shear flows — 
such as mixing layers, wakes and jets — are a class of flows characterized by the existence 
of an inflection point in the streamwise base velocity profile. For compressible flows, the 
classical inflection point is replaced by a generalized inflection point which is determined 
by the cross-space derivatives of the base velocity and density. These free shear flows 
are highly unstable. Small perturbations of suitable frequencies (or wavenumbers) in the 
flow will grow rapidly through the so-called Kelvin-Helmholtz instability mechanism. The 
instability waves will evolve into a finite amplitude disturbance which is characterized by 
streamwise-periodic regions of concentrated vorticity — - the so called Kelvin cat’s eyes. 
Although the precise mechanisms of the spreading of turbulent flows are presently still 
under study, it is generally accepted that the coherent structures play a considerable role 
in the dynamics of turbulent flows; the relative importance of the large and fine scale 
structures is still not understood all that well and the effect of the latter on the former 
has not been quantified experimentally with any degree of precision. 

Distinguishable evidence of these wave-like coherent structures in fully turbulent sub- 
sonic (unheated) round jets is convincingly displayed in the schlieren photographs of Ahuja 
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et al. (1982). They concluded that these large-scale structures, whose dimension is com- 
parable to the shear layer thickness, play an important role in controlling the dynamics 
of turbulent flows. For example, when they excited a jet by an upstream acoustic signal 
at a strouhal number, S e = F*Dj/Uf = 0.5 (F* is the excitation frequency, Dj is the jet 
diameter, and If? is the jet velocity), the jet plume was widened considerably. By using a 
photographic enhancement process, the large-scale vortices were clearly identified at some 
distance downstream. Further, they found that if the excitation frequency was doubled, 
a different instability mode was excited, the size of the coherent structure was reduced 
and the spacing between successive structures was halved. These results confirmed the 
sensitivity of the naturally existing large-scale turbulence structure to the frequency of the 
imposed disturbance. The adjacent vortical structures may amalgamate at various degrees 
of regularity as they convect downstream, thereby generating larger vortices. This pairing 
of vortices can be seen from Figure 1 of Wille’s (1963) work. When he injected smoke into 
the wall boundary layer of a circular jet nozzle, the smoke rolled up into rings of vortices 
and proceeded to pair with their neighbors farther downstream. The coalescence of the 
large scale structures was also demonstrated by Winant and Browand (1974). By injecting 
a filament of dye into a shear layer, they found that the dye rolled up into discrete two- 
dimensional (vortical) structures; farther downstream, they observed that two adjoining 
structures interact to form a single, larger vortex. Together with Laufer, Kaplan ii Chu 
(1973), they conjectured that the sequential mergings of these structures is a possible 
mechanism for the mixing and spreading of all free shear flows. 

Over the years, experiments have established that external excitations can produce 
changes in large scale and small scale motions in jet flows, although the effect of these on 
the the radiated noise field is not entirely clear. Ffowcs Williams and Kempton (1978) 
found that the increase in jet noise, due to excitation, is a direct consequence of the 
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amplified large-scale structures, while the fine-scale chaotic turbulence plays a relatively 
minor role in the noise generation mechanism. On the other hand, Morris and Tam (1977) 
argued that, although the large-scale structures are responsible for producing considerable 
changes in the fine-scale turbulence, the actual noise generation mechanism is a result of 
the latter. They provided two explanations. First, since the phase speed of an excited 
large-scale instability wave is subsonic relative to the ambient fluid, the large scale struc- 
tures are inefficient in directly generating the sound waves. Secondly, if the noise increase 
is directly caused by the large-scale structures, there would be a distinct frequency band 
centered around the frequency of an excited, most amplified, instability wave. Instead, 
jet-noise is broad-band. Nevertheless, irrespective of the mechanism through which noise 
is produced, the large-scale structures play an important part in the noise generation and 
the dynamics of turbulent shear flows. 

A distinguishing feature of coherent structures is that they are relatively permanent 
entities — in the sense that they persist for long times and distances — which co-exist with 
the seemingly completely random turbulent scales. This is in contrast to the completely 
chaotic small-scale motions which quickly lose their identities. This perseverance of the 
large scale structures offers a hope that one can manipulate the development of turbulent 
flows by modifying the evolution of the large scale structures, especially in the initial re- 
gion of jets and mixing layers. There is now considerable experimental evidence that the 
large-scale structures, and hence, the global features of turbulent flows can be organized by 
artificial excitation such as oscillating flaps, vibrating ribbons and acoustic devices (Wyg- 
nanski and Petersen 1987). By externally exciting a jet with a loud speaker upstream of 
the jet nozzle, Ahuja et al. found that not only did the radiated sound increase, but the 
strength and regularity of the coherent structures were augmented. Ahuja et al. (1982) 
presented a series of measurements for subsonic jets excited at different frequencies. They 
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have concluded that higher Mach number unheated jets require higher excitation levels 
to produce noticeable changes in both large-scale turbulence and small-scale turbulence 
intensities. One of the objectives of this research is to investigate one mechanism (i.e., a 
secondary instability theory) that might explain this experimental observation. 

The effect of excitation levels on jet mean flow velocity and turbulence intensities were 
demonstrated by Ahuja et a 1,. They found that when the excitation level was increased, 
there was a rapid decay of the mean-velocity with downstream distance. Petersen has 
shown that the mean velocity profile and the spreading rate of a turbulent, axisymmetric, 
unheated low speed jet can be altered by introducing a controlled acoustic excitation. 
This observation is important because the increase of jet spreading and the decrease of 
the mean velocity impinging on the flaps of a short takeoff and landing aircraft with an 
under-the-wing or over-the-wing externally blown flap system can lead to a reduction of 
flap loads and, therefore, a decrease in structural fatigue problems. By decreasing the 
aerodynamic and thermal loads impinging on the ground through increased mixing and 
spreading rates, the ground effects for vertical and short takeoff and landing aircraft can 
similarly be reduced. 

Stone & Mckinzie (1984) concluded that “control of flows by intentional excitation of 
natural flow instabilities involves new and largely unexplored phenomena and offers con- 
siderable potential for improving component performance.” They stressed that deliberate 
acoustic excitation — which circumvents the use of intrusive methods such as suction 
and blowing — of these coherent structures offers a promising new means of controlling 
turbulent shear flows. Most importantly, they (see also, Ahuja et al. 1986) have also 
emphasized the fact that, to date, most research work has been centered around unheated 
jet flows, usually at low speed. Therefore, future research work should be directed towards 
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the understanding and control of coherent structures in heated high speed jets. 

Control of turbulent flows via coherent structures has obvious technological significance 
for application to a considerable number of engineering devices that depend on mixing such 
as ejectors, internal mixers, combustors, jets engines and high lift devices, to name a few. 
One of the main culprits of poor ejector performance is an incomplete mixing; therefore, 
ejector performance may be improved by enhanced mixing. Enhanced mixing in internal 
mixer nozzles can increase thrust, improve fuel economy, and reduce engine noise. An 
increase in fuel/air mixing in combustors can lead not only to shorter combustors, but 
also to improved pollution control. In supersonic jet engines, the interaction of large scale 
structures with shocks is the source of shock noise (Tam & Jackson 1983); by exciting the 
jet shear layer with a high frequency, the formation of coherent structures will occur at a 
smaller scale, thereby reducing the intensity of the radiated jet noise. 

Having established the practical importance of manipulating the downstream develop- 
ment of a jet shear layer via controlling the evolution of large-scale structures, we theorize 
that the effect of manipulating these large-scale structures (or synonymously, instability 
waves) may be predicted and understood without resorting to the usual empiric si require- 
ment for the description of the fine-scale turbulence (see also, G aster, Kit ic Wygnanski 
1983). To reinforce this point, the traditional approach to turbulent flows has been known 
to be more successful in describing existing data rather than predicting new insights. 
Therefore, the goal of this research is to understand, at least in qualitative terms, how 
certain characteristics of heated and compressible jets — such as spreading rate — can be 
altered by external excitation without relying on the empirical modeling of the fine-grain 
turbulence. We will do this by a systematic examination of the stability and resonances 
of the coherent structures (assuming that they are instability waves) of subsonic heated 
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jets. 


1«2 Historical Background 

When a flow passes through a nozzle, the boundary-layer near the inlet duct at the 
jet exit quickly evolves into a free shear layer (see Figure 1) which separates the poten- 
tial core from the unperturbed ambient. When this shear layer is excited by a wave-like 
disturbance of suitable frequency, an instability wave develops in the flow. In the initial 
stages of development, where the disturbance amplitude is still small compared to the 
jet exit velocity, this wave grows in amplitude as it propagates downstream. The charac- 
teristics of this instability mode is described well by linear stability theory (Michalke 
Hermann 1982, Petersen & Samet 1988). Farther downstream, where the amplitude of 
the disturbance is of appreciable size, say a small percentage of the mean flow, nonlinear 
effects become important, and the wave (or coherent structure) reaches a finite-amplitude 
saturation. An important objective of this work to understand and to promote the (reso- 
nance) conditions under which this instability wave can destabilize a subharmonic mode. 
Physically, this is important since a subharmonic disturbance can generate vortex pairing; 
this pairing is, at least partially, responsible for jet spreading and mixing. 

The existence of large-scale instability waves has now been acknowledged by an over- 
whelming number of observations; an extensive review of coherent structures in excited 
shear flows is given by Wygnanski & Petersen (1987). Crow and Champagne (1971) pro- 
vided the first comprehensive work on the response of a circular incompressible unheated 
jet to a controlled axisymmetric excitation of a certain amplitude and frequency. By 
measuring the velocity fluctuations along the jet axis, they found a “preferred” mode at 
which an excited instability wave reached a maximum amplitude. This preferred mode 
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has a frequency of f p = O.ZUf /D* where Uf is the jet exit speed, and D * is the diameter. 
They concluded that this mode is an instability wave most capable of reaching a large 
amplitude before saturating. The significance of this preferred mode lies in the fact that 
it is the most dominant and frequently occurring of all large scale structures at the end 
of the potential core. Michalke (1971) determined that the measured phase velocity of 
the coherent structure, found by Crow &z Champagne, agreed reasonably well with that 
obtained from a linear stability theory based on a turbulent mean flow. By using phase 
locked measurements of controlled excitations in an axisymmetric cold jet, Petersen ic 
Samet (1988) concluded that the preferred mode is, in fact, a Rayleigh instability mode, 
provided that the stability analysis is based on the local shear layer thickness and the 
measured mean velocity profiles. As a result of this finding, they concluded that the 
preferred mode in, both excited and naturally occurring, jets is a shear layer instability 
mode rather than a global instability mode involving the entire jet column. This lends 
considerable credibility to the use of stability theory to understand the nature and the 
evolution of the quasi-deterministic large scale motions associated with the appearance of 
coherent structures in both unexcited and excited jets (see also Gaster, Kit it Wygnanski 
1985). 

There have recently been a number of theoretical works on the inviscid stability charac- 
teristics of axisymmetric jets, although studies of either compressible or heated jets remain 
scarce . Batchelor <k Gill (1962) pioneered a theoretical analysis on the inviscid stability 
characteristics of a top hat velocity profile which characterizes the mean flow close to the 
jet exit. They found that the top hat velocity profile is unstable to a small disturbance for 
all axial and azimuthal wavenumbers. Their analysis is supported by Plaschko (1979) and 
Cohen (1986). Stability theory has now been successfully extended in various ways and 
applied for jets at near ambient temperatures and at low jet Mach numbers (for example, 
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Cohen 1986). 

The stability characteristics for hot jets, which are of technological interest, have re- 
ceived relatively little attention. Michalke (1971) and Michalke <k Hermann (1982) are 
among the few who investigated the stability characteristics of subsonic heated jets by 
the linear, inviscid and parallel flow stability theory. They found that the local growth 
rate of spatially growing waves increases as the total temperature of the jet is increased. 
This means that for heated jets there is a more unstable coherent structure, resulting in 
the widening of the mean jet velocity profile. Recently, Sohn (1986) found that cold jets 
are convectively unstable, whereas hot jets with an exit temperature 1.5 times the ambi- 
ent temperature are absolutely unstable. By using a “quasi-linear* model, Ahuja et. al. 
(1982) concluded that the small-scale turbulence, which is induced by the passage of ex- 
cited large-scale structures, is responsible for the jet spreading and mixing for both heated 
and unheated jets. Their theoretical and experimental studies have shown that both jet 
noise amplification and jet mixing decrease as the Mach number increases. The effect 
of heating is predicted to increase the turbulence intensity and to widen the jet width. 
They pointed out that their theoretical results for the heating effect on jets are inconsistent 
with limited experimental data and, therefore, are inconclusive. The “quasi-linear model,” 
which requires empirical constants and certain turbulence closure assumptions, provides 
little physics governing the role of coherent structures in turbulent jet flows. Hence, the 
mystery and intricacies of coherent structures in heated jets have yet to be systematically 
unravelled. 

Over the past few years, it has been established that linear inviscid stability theory can 
predict with surprising accuracy the phase velocity and the passage frequency of coher- 
ent structures, as well as the lateral distribution of the perturbations in excited jets (e.g., 
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Crighton k Gaster 1976, Plaschko 1979, Strange k Crighton 1983). The convincing agree- 
ment obtained by Gaster, Kit k Wygnanski (1985) between experimental data and the 
linear inviscid stability theory for slowly diverging shear layers reconfirmed the preceding 
remark. Consequently, coherent structures can be viewed as (inviscid) instability waves 
of the mean flow, and fine-grain turbulence plays a passive role in the development of the 
coherent structures. They stressed that from the theoretical point of view, the large-scale 
structures offered a hope that certain characteristics — such as mixing enhancement and 
suppression — may be predicted and understood without resorting to the empirical re- 
quirement for the description of the fine- grain turbulence. 

Although the linear parallel flow instability theory can describe accurately the local 
stability characteristics of the initial stages of the development of a large-scale instability 
wave, it is unable to predict the overall streamwise amplification of the disturbance. By 
retaining weakly nonparallel terms and using the method of multiple scales, Crighton k 
Gaster (1976) gained substantial improvement in the total streamwise growth, although 
there still exists discrepancies between their numerical results and experimental data. Pe- 
tersen k Wygnanski (1987) attributed these differences to the neglect of the nonlinear 
terms in a linear analysis. 

Another deficiency of linear (parallel flow) theory is that it is unable to predict the 
growth rate of a subharmonic disturbance in a periodic base flow. In the case of a planar 
mixing layer, this growth rate was first calculated by Kelly (1967) through the so-called 
\ineai secondary stability analysis. His temporal analysis predicted that a periodic funda- 
mental component superimposed on a parallel mean flow destabilizes a subharmonic dis- 
turbance whose phase velocity matches that of the fundamental. Purely numerical work 
on the instability of a periodic flow has been carried out by Pierrehumbert and Widnall 
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(1982) who studied the linear (temporal) secondary instability of a row of Stuart vortices. 
Their results are in reasonable agreement with Kelly’s analysis. Cohen (1986) extended 
Kelly’s analysis to spatially growing waves in a planar mixing layer. While Kelly decoupled 
the predominant fundamental wave from the growth of the subharmonic, Cohen consid- 
ered the exchange of energy between a fundamental disturbance and its subharmonic as 
both waves propagate downstream. Additionally, Cohen established subharmonic and az- 
imuthal resonance conditions in low speed and unheated jets. In particular, he has shown 
that the nonlinear interaction between two azimuthal modes can produce a third mode 
which is initially absent in the flow. His theoretical results are in close agreement with 
his experimental studies. More recently, Monkewitz (1988) generalized Kelly’s analysis 
and examined vortex pairing in more detail in incompressible mixing layers. By incorpo- 
rating explicitly the weakly nonlinear interaction terms between two instability modes in 
his method of multiple scales analysis, he found that the development of the subharmonic 
depends crucially on its phase relation with the fundamental. He has also shown that a 
critical threshold fundamental amplitude is required for the subharmonic to become phase 
locked with the fundamental; this presence of the fundamental leads to an enhancement 
of a subharmonic growth rate. 


It is perhaps appropriate at this point to mention briefly that the linear secondary 
stability mechanism is also important in the studies of the laminar-turbulent transition 
in walled-bounded boundary layer flows. By investigating the temporal instability of 3-D 
disturbances in a spatially periodic base flow consisting of a Blasius velocity profile and a fi- 
nite amplitude 2-D Tollmien-Schlichting wave, Herbert (1984) found that the subharmonic 
resonance mechanism is the principal route to transition in low disturbance environments; 
hence the designation H-type breakdown. His numerical results agree closely with the 
experimental data of Saric (1983). 
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1.3 Scope of Present Research 

The purpose of this research is to examine in detail the stabilities and resonances 
of coherent structures — assuming that they are synonymous with instability modes — 
in heated subsonic round jets. We choose temporal stability theory in favor of the spatial 
theory because the validity of a spatial theory was unclear at the initiation of this work (see 
remarks on convective instabilities of heated jets in Chapter 2). We believe that a temporal 
analysis of heated jets will contain all the essential physics and qualitative trends, although 
it is certainly true that the corresponding spatial analysis may provide a better agreement 
with experimental data. Of course, spatial instability arises naturally in convectively 
unstable flows when periodic excitation is present. The point is that the dispersion relation 
(i.e., the eigenvalue relation) is an analytic function of the complex wavenumber, and 
spatial instability is merely the analytic continuation of temporal instability. In addition, 
there are a number of conceptual difficulties associated with spatial instability — such as 
the calculation of the mean flow distortion — which we wish to circumvent through the 
use of temporal analysis. The present investigation is divided into the following two tasks: 

1. Although the linear parallel flow instability theory of incompressible flows is presently 
understood and documented, the stability analysis of compressible and heated round 
jets is still rather sketchy. In the first part of this work, we focus on the systematic 
determination of the individual role of Mach number and temperature ratio in affect- 
ing jet instabilities. A good understanding of the primary instability characteristics 
is also a prelude to secondary instability calculations. In this work, an analytic mean 
flow (tanh shape) which provides a good local approximation to that obtained from 
measured experimental data is used in our stability calculations. 

2. The principal task of this work is to study the linear subharmonic instability of a peri- 
odic flow consisting of the steady mean flow and a small but finite amplitude primary 
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wave. In order to make the physics governing this instability more transparent, we 
first establish the conditions — the so-called resonance conditions — which provide 
an effective mechanism for the destabilization of an instability wave. The resonant 
interactions between two waves allow us to determine the relevant parameters to be 
used as input for the secondary instability calculations. 

We carry out the “secondary” instability analysis using the methods of multiple scales 
and a generalization of the normal mode analysis on a streamwise periodic flow. It 
may be noted that the secondary instability of round jets is, to date, still new, in 
the presence of heating and compressibility. The growth rate of the subharmonic 
disturbance as a function of the amplitude of the primary wave for different jet 
parameters — such as Mach number, temperature ratio, shear layer thickness and 
modal numbers — are determined. 

In summary, this research will provide a deeper understanding of the relevant param- 
eters (e.g., wavenumber of excitation, modal number, Mach number and jet temperature) 
that affect the secondary instability and the physical mechanisms governing vortex merg- 
ings in subsonic heated round jets. 
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Chapter 2 

GOVERNING EQUATIONS AND PRIMARY 
STABILITY PROBLEM 


2.1 Governing Equations 


The governing equations for a viscous, heat conducting and compressible ideal gas 
in dimensional form (characterized by a superscript asterisk) are the conservation of mass, 
momentum, energy, and the thermodynamic equation of state. They are 

§£ + V*.(pV) = 0 (2.1) 

W = “ V V + + V* • r’ (2.2) 

P * C * V 1FF = ~ P * V * ' v * + V * ' (**VT*) + & (2.3) 

p * = p’STT' (2.4) 

where p*, v*, p* and T* represent the density, velocity, static pressure and the temperature, 
respectively. A point in space is denoted by the Cartesian coordinate x*, V* = <9/ 3x* 
and t* stands for time. The convective derivative is defined by 


D* 

D't * 



V*, 


(2.5) 


g* is the constant gravitational acceleration, C* is the specific heat capacity at constant 
volume and is the gas constant. In (2.3), Fourier’s Law of heat conduction is assumed; 
the heat flux is given by (— k* V + T*). These equations may be found in Bird, Stewart and 
Lightfoot (1960). 
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Dissipative processes in the fluid are due to the viscosity, p*, and the thermal con- 
ductivity k* . For a Newtonian fluid, the shear stress tensor, r* , and the “dissipation 
function”, $*, are defined by 

r* = 2/i*S*-H A1 *(V*-v*)I (2.6) 

= r * : S’ = 2/i*S* : S* - -v*) 2 (2.7) 

3 

where 

s* = 1[VV + (V*v*) T ] (2.8) 

M 

S* is the rate of strain tensor, I is the idemtensor and (*) T denotes the transpose of a 
tensor or matrix. 


2.2 Nondimensionalizing the Governing Equations and 

Simplifying Assumptions 

Our principal interest is in the study of certain types of instabilities associated with 
compressible and heated round jets. For this purpose, the governing equations (2.1) - (2.4) 
can be made dimensionless by using the centerline values at the jet exit of the base flow 
(also called the mean flow or the unperturbed flow). We nondimension alize all length, 
velocity and time scales by i?*, Uj and R*/Uf y respectively — where R* is the radius of 
the jet. The dimensionless physical variables become: 



31 


and with the above normalizations, the governing equations (2.1) - (2.4) become 

^ + V-pv = 0 (2.9) 

v, + £*.. + (210) 

= -ft - 1 ) 7 M/pV • v+ j^-V-kVT + 7,7 ~ ' )M) » (2.11) 

1 M]p = pT ( 2 . 12 ) 

where 

S=i[Vv + (Vv) r ] (2.13) 

r = 2 M S-|^(V.v)I (2.14) 

$ = r:S (2.15) 


and k is the unit vector in the direction anti-parallel to the gravity. The Reynolds number, 
Re } Mach number, My, Froude number, Pr, and Prandtl number, Pr, are given by 


Re = 




uf 

Fr — -V. 

g'R' 


Pr = 


_ Mjc; 


where is the unperturbed sound speed at the jet exit, C* is a constant representing the 
specific heat at constant pressure, 7 = constant = C*/C*, and g* is a constant representing 
the magnitude of the gravitational acceleration. The nondimensional convective derivative 
is given by 

£ = !+v.V (276) 


In order to make any progress with the relevant equations, a large number of assump- 
tions must be made which are discussed in books on stability theory (Lin 1955, Betchov 
and Criminale 1967, Drazin and Reid 1981). To begin with, the last term in ( 2 . 10 ) is 
ignored because gravitational effects are unimportant at large Froude numbers; this is a 
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standard assumption in aerodynamics, for example, when R* = 5 cm, U* — 5000 cm/sec, 
p* = 1.3 x 10" 3 gm/ cm 3 , /ij = 1.8 X 10" 4 dynes sec/cm 2 (air at normal conditions), Fr = 
5 x 10 s , Re — 1.4 x 10 5 , Pr = 0(1). In instability theory, it is convenient to decompose 
the dependent variables into two parts: the base flow and the perturbations. At high 
Reynolds numbers, the base flow (i.e., unperturbed flow) will satisfy the boundary layer 
equations so that the unperturbed pressure will be a constant (approximately) throughout 
the flow field, and the jet will “diverge” with downstream distance on a long length scale 
of (i?*/?e5). Since the characteristic wavelength of the instability wave is on the order 
of R* y the divergence of the base flow may be ignored for the purpose of calculating the 
instability wave. In other words, the base flow is assumed to be locally parallel. The 
precise form of this flow will be specified in the next section. 

At high Reynolds number (say, above 10 3 ), the effects of viscosity and heat conduc- 
tivity are negligible for the instability wave, except possibly for neutral modes find their 
derivatives in the cross-space. For this reason, we shall retain a highly simplified version 
of the viscous and conduction terms: these terms are only important in a thin layer of 
thickness h = O (iJe -1 / 3 ), centered around the generalized inflection point. In this layer, 
the so-called the critical layer, p and k are constants to a high degree of approximation 
so that V • r may be replaced by /i c V 2 v and V.(fcVT) by k c V 2 T y where the subscript 
c denotes the value of a quantity in the critical layer. As a result of this approximation, 
the effects of heat conductivity and viscosity are correctly represented only in the critical 
layer where they may be important. Outside the critical layer, these effects are incorrectly 
represented in our analysis, but this is immaterial because they are unimportant. 

Further, when the disturbances are small, the largest terms in $, in (2.11), are propor- 
tional to the product of the velocity gradient [which is O(l)] and the perturbation velocity 
gradient [which is at most O ( 8/h), where 6 <C 1 is the characteristic magnitude of the 
disturbance]. On the other hand, the perturbation temperature gradients in the conduc- 
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tion term of the energy equation scale as S/h 2 because of the presence of the Laplacian 
and are, therefore, much more important in the critical layer than Consequently, the 
viscous dissipation, will be ignored here. 

In summary, viscosity and heat conduction are included, but only to the extent needed 
to eliminate any possible singularity in the instability modes arising from the reduced 
inviscid and nonconducting equations. Although the viscosity and conductivity of a fluid 
vary with temperature, these variations will be negligible across the thin critical layer 
where they play an appreciable diffusive role. Consequently, the viscosity and conductivity 
are assumed to be uniform throughout the fluid. It may be remarked that outside the 
critical layer, the viscous and the conduction terms are negligible relative to the convection 
terms. With the preceding assumptions, the governing equations are approximated by 
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2.3 The Base flow and Perturbation Equations 

We shall now formally separate the disturbance from the base flow by writing 


P^> Po + P 


v — ► v 0 + v 


( 2 . 20 ) 


p-*Po + P 

where the subscript 0 hereafter refers to the base flow, and p, v, p are perturbation 
quantities . After substituting (2.20) into (2.17) - (2.19) and collecting terms, we obtain 
a set of equations for the base flow and for the perturbations. The latter set is given 
vectorially as 


P o- 


Dq P 
Dt 

D 0 v 


+ p 0 V ■ v + v • Vp 0 = -V • (pv) 


Dt 


+ pqv ■ Vv 0 + Vp - -j^V 2 v = -p-^7 


Dv 

Dt 


^ + 7 PoV • v - -=^=-V 2 (— - = -v • Vp - 7 pV • v 

Dt PrRe po p 0 

+ ^L_v 2 (^-^ + .-) 

PrRe Po 


( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 


where, 


Do _ d , d 
Dt dt 00 dx' 


(2.24) 


and the dots stand for terms for the product of perturbations that are cubic and higher. 
In deriving (2.21) - (2.23), we have made the assumption that the base flow is parallel; 
this implies that 


Vo = [uo(y,2),0,0] 

(2.25) 

Po = My,z),o,o)] 

(2.26) 

po = constant 

(2.27) 
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where x = (x,y, z) is a cartesian coordinate system with x pointing dong the axis of 
the jet. Although the numericd calculations and the theoretical developments are carried 
out in a cylindrical coordinate system (x,r, <)>) with y = rcos^, z = rsin<£, it is more 
compact to write down the perturbation equations in a Cartesian system. Note that Vpo = 
(0, dpo/dy , dpo/dz) so that v ■ Vpo does not involve the x-component of the perturbation 
velocity. Similar remarks hold true for other terms of this type. In order to generate 
specific results, the mean flow is taken as 


where 


and, 


«o = uo(»0 = i{l - tanh[(r - -)/i 0]} 

At T 


Po = Po(r) = 


n 

1 + (T. - l)uo + ^M/r.uo(l - tio) 
1 

po = constant — — 

1 Mf 


fOO 

9 is the dimensionless momentum thickness = / uo(l — uo)^ 

Jo 

p— . . . r 0 (o) 

is the temperature ratio = ■ / \ 

F To(oo) 


To is the base temperature = To(y,x) = To(r) 


(2.28) 

(2.29) 

(2.30) 


Note that the steady unidirectiond velocity profile (2.28) is not an exact solution of the 
Navier-Stokes equation. It is, however, a good local approximation to that obtdned from 
the measured experimentd data. This mean velocity profile, which has dso been used by 
other investigators (e.g., Michalke & Hermann 1982), is capable of representing the top- 
hat profile at the jet nozzle exit and the fully developed Gaussian profile at around the 
end of the potentid core as well as the profiles between these two extremes. This can be 
done by varying the parameter 9 (momentum thickness) from (near) zero to a vdue which 
is a substantid fraction of unity (i.e., the nondimensional jet radius). The parameter 9 
can thus be used to characterize different velocity profiles at different axial locations. 
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For compressible ideal gas flows, the mean density profile, (2.29), is obtained from 
the mean velocity profile via the Busemann-Crocco law, which is valid for high Reynolds 
number flows and when the Prandtl number is unity (Schlichting 1979, p. 330). For a 
given velocity profile, uo, the mean density distribution is governed by three dimensionless 
parameters: Mach number My, temperature ratio T,, and the isentropic exponent 7 . In 
this work, the isentropic exponent 7 = 1.4, and T* is taken to be greater than one, which 
means that the jet is hotter than the quiescent environment. For example, when T, = 2, 
the jet is twice as hot as the environment. 

Figure (2.1a) shows the dependence of velocity profiles on 9 while the effect of 9, Mj, 
and T* on the mean density distributions are shown in Figure (2.1b). It is interesting 
to note that in Figure (2.1b), T, has a much greater influence than the (subsonic) Mach 
number My upon the base density profile, and therefore, is expected to produce a greater 
effect on the instability. 


2.4 Primary Linear Stability Analysis 

Although our principal goal is the study of the secondary instabilities and sub- 
harmonic resonances in a compressible and heated jet, we begin our discussion with the 
classical Unear instabilities of a jet. A knowledge of these linear instabilities is needed in 
order to study the secondary instabilities and resonances. Perhaps it is worthwhile, first, 
to say a few words about the parallel flow assumption for the mean flow in view of some 
relatively recent developments. 

Petersen & Samet (1988) have found that the local stability characteristics (e.g, eigen- 
functions, phase velocity , Reynold stress, etc.) of an instability mode developing on a 
jet column can be predicted accurately from linear viscous stability theory based on the 
measured mean velocity profiles, which are again assumed to be locally parallel. They 



37 


found that the preceding observation is valid even when the local streamwise velocity dis- 
turbance reaches a level as high as 24% of the jet speed, when the flow is unquestionably 
nonlinear. The last remark lends considerable credibility to the generally accepted claim 
that even in “nonlinear” flows, the mode shape of the fundamental is given quite accu- 
rately by the linear mode. A large body of literature has been built on this assumption, 
beginning with the pioneering work of Ko, Kubota and Lees (1970). 

In fact, this finding of Petersen <k Samet is important in the instability analysis of 
a jet which has a finite- amplitude wave- like structure developed on its column. This 
structure is approximately both spatially and temporally periodic. In order to represent 
the fundamental component of this periodic flow, we shall use a linear instability mode. 
This makes the study of secondary instabilities much easier because we do not have to 
compute a finite amplitude primary disturbance. In our secondary instability analysis, 
this amplitude is a given quantity; the mode shape is given by linear theory. 

Unfortunately, a shortcoming of linear instability based on parallel flow assumption 
is that it predicts poorly the total streamwise amplification of an instability mode. By 
retaining weakly nonparallel terms in an axisymmetric round jet, Crighton ic Gaster (1976) 
[see also Gaster, Kit and Wygnanski 1985] gained a significant improvement on the total 
growth rate, although there still exists some unsatisfactory discrepancies between their 
results and experimental measurements. According to Wygnanski As Petersen (1987), 
these discrepancies are attributed to the neglect of the nonlinear interaction terms in the 
linear parallel flow analysis of a single wave train. Petersen et al, have, in fact, dispelled 
the skepticism surrounding the predictive capability of a parallel flow stability analysis 
for the mode shape. While there is no doubt that an accurate determination of the total 
growth rate still requires further studies on nonlinear analysis, it is not a concern in this 
work. The quantities of interest here are the local values of the phase velocity, growth 
rate and the shape of a specific mode, and that can be described well by linear parallel 
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flow analysis. 

Because of the effects of compressibility and heating in a cylindrical geometry, the 
algebraic aspects of the analysis in this thesis are very lengthy. In order not to burden 
the reader, we use a compact operator notation; the specific elements of these operators 
may be found in the various appendices. 

To recall, we use a cylindrical polar coordinate (x, r, ^), where x is along the axis of the 
jet and (r, <f>) are the radial and azimuthal coordinates; the components of the perturbation 
velocity are v = (u,v,u;) in this cylindrical system. To represent the entire perturbation 
field, we use the five-dimensional vector 

F = S f 4 = S (u,v,w,p,p)l 


After substituting the above equation into the (nonlinear) perturbation equations, (2.21) - 
(2.23), and neglecting all quadratic and higher order terms in 8 , we arrive at the linearized 
disturbance equations 

A b f b = 0 (2.31) 

Here A* is a linear operator whose action on is defined in detail in Appendix A. The 
various terms in A* are functions of the base flow alone; A* has variable coefficients which 
depend only on the radial coordinates. Here 8 is a small parameter which measures the 
magnitude of a typical disturbance with respect to the base flow. Since A b depends on 
r alone, we may extract the dependence of the disturbance on the other coordinates via 
Fourier (space) and Laplace (time) transforms. Because of the linearity, we consider a 
“single mode” (or wave component) and write 


f b = f(r)e*( ax+ ^e ff * + complex conjugate; * — \Z~"l 
f (r) = [u(r),t;(r) ) «;(r),p(r),p(r)] r 


(2.32) 


where a and j3 are the wavenumbers in the axial and azimuthal directions, respectively, and 
a is the complex growth rate. The disturbance is completely described by a, /?, a and the 
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complex amplitude function (also known as the mode shape) f(r). The modal form (2.32) 
is simply a single mode of a Fourier-Laplace component representing the wavenumber and 
frequency component of an arbitrary disturbance. This means that each individual mode 
evolves independently in the course of time, according to the the linearized equations. 
Because of the periodicity in <£, f3 must be an integer (0, 1,2 ... etc), and it is also called 
the mode number. In temporal theory, a is real, and a is complex. The real part of a 
determines the temporal growth rate of a disturbance. The imaginary value of a is the 
radian frequency. In spatial stability, a is usually written as — iu where oj is the given 
real frequency, but a is complex. The spatial growth rate of a disturbance is given by 
the imaginary value of a. An approximate relation between the growth rates of these 
two forms of instability can be obtained by Gaster’s transformation (Gaster 1962) under 
certain restrictive assumptions — such as for small growth rates. 

Although there still exists some questions on the application of temporal or spatial 
theory (especially for nonlinear flows), it is now generally accepted that temporal analysis 
is meaningful for absolutely unstable flows while spatial analysis has physical meaning 
for convectively unstable flows. The basic concepts of absolute and convective instability 
are discussed in detail by Huerre As Monkewitz (1985). Suppose a jet is perturbed by 
a disturbance which is localized in space and is impulsive in time. Some time after the 
triggering of this disturbance, a wavetrain of finite extent evolves on the jet. If both the 
leading and trailing edges of this wavetrain propagate downstream , the flow is said to 
be convectively unstable. Otherwise, the flow is absolutely unstable. In an absolutely 
unstable flow, the trailing edge of the wavetrain propagates upstream and the group 
velocity is negative in some region of the wave. Because of this, disturbances which are 
generated downstream of the nozzle in an absolutely unstable flow can actually interact, 
through instability waves, with the nozzle to produce “feedback” or “resonance’ loops. 


Temporal stability is also applicable to flows which are required to be periodic in space. 
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Such flows may not actually exist in practice, but are at least fairly close approximations 
to those that do. Once a flow is forced to be periodic in space, its natural spatial evolution 
is inhibited and the (linear) instability will appear as an exponential growth in time. For 
this restricted class of flows, temporal instability is always valid. Since at the beginning of 
this study, the concepts of convective and absolute instabilities were not clarified, we opted 
for a study which is based on temporal instabilities and spatial periodicity. In retrospect, 
this was a good choice because it is now known that hot jets with T* >1.5 may be 
absolutely unstable (Sohn 1986); for such flows a study based on spatial instability would 
be meaningless. 

The modal representation (2.32) of a disturbance transforms the system of linearized 
partial differential equations, (2.31), into a system of ordinary differential equations (here- 
after ODE’s) in r. These equations can be conveniently put into the matrix form 

Aif = aBif (2.33) 

where the elements, which depend on the wavenumbers (a,£), of Ai and Bi are given in 
Appendix B. The principal task is to solve (2.33) for a given (a,/?) and suitable boundary 
conditions in order to obtain the eigenvalue from the dispersion relation a = <r(a,/?). We 
now discuss the boundary conditions associated with (2.33). 


2.5 Kinematic Boundary Conditions on the Linearized 

Disturbances 

The appropriate boundary conditions for a disturbance superimposed on a subsonic 
jet require that f remains finite on the jet axis and vanishes at infinity. They are, 

f — ► 0 as r — ► oo (2.34) 


f remains finite at r = 0 


(2.35) 
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In fact, some of the boundary conditions to be satisfied at the origin are purely “kinematic” 
in nature, and therefore, are independent of the viscosity (Batchelor and Gill 1962). In 
order for a disturbance to remain single- value at the origin, these kinematic boundary 
conditions yield 

u(0) — p(0) = p(0) =0 for /? ^ 0 (2.36) 


and, 


v(0) = w(0) = 0 for $ ^ I 


(2.37) 


Other boundary conditions at the origin, which cannot be determined from kinematic 
consideration alone, are derived from their asymptotic behavior developed in the next 
section; note that not all the boundary conditions are given by kinematic considerations. 


2.6 Asymptotic Behavior of a Disturbance on a Uniform 

Base Flow 

The numerical solutions of hydrodynamic stability problems in an unbounded do- 
main, such as in round jets, are often solved over a finite domain. Consequently, the 
analytical structure of a disturbance for large r and small r is required. Because the base 
flow, uq and /?o> approaches a constant value, the asymptotic behavior can be obtained 


explicitly by solving a series of Bessel-like equations. In fact, for r — ► 0, we can write 

u = ciu p + c 2 u v (2.38) 

t; = civ p + c 2 v v (2.39) 

w ~ ciWp + c 2 w v (2 40) 

P = cipp (2.41) 

(2.42) 


P = ci p p 
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where 

u p = y-/fl(A p r) 

A p 

V p = Ip{ r ^p) 

Wp = ^b( rX p) 

[A* -a 2 . , 

Pp =|-^--p°.a(u°-C ) |^— 

a 2 

Pp = Po(~r ~ A p )i/j(A p r) 

A p 

^ _ a[l - ^ 0 («o - C) 2 M 2 + ia(«o - C)Mj/Re]i 
P [1 + ia(u 0 - C)Mj/Re]? 

u v = ^-//j(A„r) 

v v = I'p{\ v r) 
iff 

w v = —I t j(A„r) 

A u = [a 2 + iapo(u 0 - C)Re]? 

where ci and c 2 are constants, the primes denote differentiation with respect to the argu- 
ment in parenthesis, and the subscripts p and v refer to the so-called “pressure mode” and 
“viscous mode” , respectively. Here Ip is the modified Bessel function whose order is the 
azimuthal wavenumber, ft = integer, and C = itr/a is the complex phase speed. Clearly, 
the kinematic boundary conditions are contained in the above asymptotic expressions. In 
order to maintain a highly accurate finite difference numerical procedure, the asymptotic 
boundary conditions are enforced. This not only fulfills the requirement of the value of a 
specific quantity at a specific point, but also the manner it behaves (e.g., in powers of r). 

Similarly, as r — ► 00, the same asymptotic results hold provided that Ip is replaced 
by the modified Bessel function of the second kind, Kp (Abramowitz <k Stegun 1972), 
and uq and po are replaced by their values at r = 00. Note that in all likelihood, the 
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viscous solution is unimportant for large r because the argument of Kp is large (i.e., 
A „r ~ (aRe)%r » 1). 

2.7 The Behavior of Disturbances in the Critical Layer 

Although we have retained the most important effects of viscosity and heat conduc- 
tivity in the critical layer, it is instructive to examine the nature of the flow in this layer 
when (formally) Re = oo in equation (2.33). It is known that under certain conditions 
— to be made more precise momentarily — the instability modes will become singular at 
the critical point. It is desirable to know the “strength” of this singularity for numerical 
purposes, and to see under what conditions at least some of this singularity may be elim- 
inated, even when Re = oo. The latter condition gives us a generalization of the Rayleigh 
inflection point criterion to compressible and heated flows in round jets. We emphasize 
that, strictly speaking, there will be rjo singularity in the modes because the “diffusive” 
effects of viscosity and conductivity have been retained. However, at large Re, we will 
have a “near” singularity that does affect the accuracy of the numerical solutions. 

By substituting Re = oo and a — -iaC in equations (2.33), the linearized, inviscid 


and nonconducting disturbance equations can be obtained as 

p 0 ia(u 0 - <?)u + Pou f 0 v + tap = 0 (2.43) 

#}ia(u 0 - C)v + p 1 = 0 (2.44) 

£o*<*(tio - C)w + = 0 (2 45) 

tau + ~ + v’ + + tar(uo - C)p = 0 (2 46) 

«a(u 0 “ C)p + Po v + Po(~ + v 9 + + tau) = 0 (2 47) 


where primes denote differentiation with respect to r. After some algebraic manipulations 
in favor of the pressure, the above equations can be reduced to a second order ODE with 
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variable coefficients (essentially a Rayleigh equation) of the form 

+ (i - -&B - + [•*(«. - c)VoM= - £ - »■]„ = o 

dr 2 r uq — C po dr J r* 


(2.48) 


When the instability wave is neutral, the differential equation (2.48) has a regular singular 
point at the point r = r c where uo(r c ) = C (C is real for a neutral wave); the subscript c 
hereafter refers to the value of a variable at the critical point. In the neighborhood of r c , 
the coefficients of the differential equation can be expanded in terms of powers of r — r c . 
For example, 

u 0 = c + u' c (r - r c ) + u"(r - r c ) 2 /2+ , 

(2.49) 

po = Pc+p' e (r-r c ) + p"(r-r c ) 2 / 2 H 

where u' e = Uo(r c ) etc. Substituting (2.49) into the differential equation (2.48), we obtain 
the approximate differential equation around the critical point for p, i.e., 

2 


p" - [ h Di + Di{r — r c ) H ]p' + [E\ + Ei{r - r c ) H ]p = 0 


r — r, 


(2.50) 


where 


and 


D = Pc< + p'c< _ I 
1 Pc u'e r c 

_ 2u” K".P?_L 

2 ~ 3< + 2«? p c + f l r? 


El = -(«’+£) 

r c 


r c 

According to the method of Frobenius (Bender ic Orszag 1978), equation (2.50) has two 
linearly independent solutions of the form r = r c , 


Pa = (r~ r c ) 3 ffi(r ~ r c ) 


(2.51) 


Pb = 9i{r- r c ) + Kp a log(r - r c ) 


(2.52) 
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where 

9i{r ~r e ) = l + <*i(r - r e ) + a 2 (r - r c ) 2 4 

gi{r - r c ) = 1 + bi(r - r e ) + b 2 {r - r c ) 2 4 

The coefficients a's, 6's and K are found by substituting the Frobenius series (2.51) - (2.52) 
into the differential equation (2.50) and equating like powers of (r — r c ): 

ZD! 

°i = j 


a l = 


ZD 2 + 4ajDi — Ei 

10 


6i =0 

h - El 

b2 ~ T 


K = 


D\E\ — £?2 


The gener 2 il solution of the differential equation (2.50) is of the form 


P = Ai p a + \2Pb 


(2.53) 


where Ai and A 2 are constants. The solution p a is analytic, while is a multi-valued 
function because of the logarithmic singularity unless K = 0. For the particular case of 
an unheated jet where po is a constant, the condition K = 0 reduces to 


d . ru 

~1Z \ 


0 \ _ 


) = 0 at r — r c . 


(2.54) 


drV 2 + « J »- 2 ' 

Condition (2.54), which was established by Batchelor A: Gill (1962) to be a necessary con- 
dition for the existence of an simplified inflectional type of disturbance in an incompressible 
unheated round jet, is, in fact, valid even in the presence of compressibility. Furthermore, 
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the familiar generalized inflection point theorem for the existence of instability modes in 
planar mixing layers, derived by Lees & Lin (1946), i.e., 


±(pou'o)=0, (2.55) 

at y = y c where y is the transverse coordinate, can be recovered by letting l/r c — ► 0 in 
the condition K = 0. 

The behavior of the velocity and density perturbations in the vicinity of the critical 
point is obtained by substituting (2.51) and (2.52) into equations (2.43) - (2.47). The 
corresponding two independent solutions near the critical point, r = r c , are 


3i , 

Va = (r - r e ) + 

u'cPc* 


3 pc*" + (6p' c - 8ai p e )u' c 


2aip 2 c u 


(r ^ r c ) 2 + 


2t &2 


v 4 = H 1- Xv a log(r - r c ) 

apcW c 


«« = — 


3 , (V. ~ 4fl lPc) 


+ 


a 2 p c u' e 1 a 2 />? u ' c 


~ r e) + 


U t = 


r|c»! 2 p (: u' (r — r c ) 


H 1- Ku a log(r - r c ) 


(2.56) 


tv a — -■ 


0 


otr c p c u' c 


j{r-r e f + 


m = - 


0 


u' c Pc r M r - r c ) 


H b A'tUalogfr - r c ) 


_ &r c u'J - (4air e + 3)u' , 

Pa — o ig I 


a 2 r e u' e 3 


J W ± A.K + jffi? - feH + *)K | + ... + ^ Mr -„) 

( <* 2 r*u*(r - r c ) J 

The above equations (2.56) show that the logarithmic term is present in all perturbations 


unless the generalized inflection point criterion is satisfied; that is K = 0. In addition, the 
leading order term in the pressure fluctuation and velocity perturbations of compressible 
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flows has the same analytical behavior as in incompressible flow. For example, the leading 
term for u and w has a simple pole unless the azimuthal mode number is zero, irrespective 
of compressibility. The density fluctuation, or the temperature fluctuation, which is a new 
quantity due to compressibility, has generally a simple pole at leading order. Consequently, 
for a neutral mode, the p and v perturbations are finite, the u and w modes becomes 
infinite for ^ 0, and the p perturbation is generally infinite. It is precisely the singular 
nature of the perturbations at the critical point that creates difficulties in the numerical 
calculations for neutral modes. With the inclusion of dissipative effects, some of this 
difficulty disappears, although a very fine grid is needed in the vicinity of the critical layer 
to resolve the fluctuations. 


2.8 Summary 

In summary, the ODE’s, which is symbolically expressed by (2.33), is solved by two 
independent numerical schemes (to be discussed later) with suitable boundary conditions 
on the jet axis and in the quiescent medium at infinity. Nontrivial solutions exist only for. 
a certain set (a, ,5, a), or equivalently in the more usual notation, 

ct-<t{ol,P) (2.57) 

The above equation is called the dispersion relation whose real part (if positive) indicates 
instability and whose imaginary part is the frequency of the wave. 



2.6 
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49 


Chapter 3 

FORMULATION OF THE SECONDARY INSTABILITY 

PROBLEM 


We begin our discussion in this section with a physical description of how an 
instability wave evolves on a jet and how the flow field may be approximated at a couple 
of diameters downstream of the jet nozzle. 

When a jet is perturbed by an external excitation of a suitable wavenumber or fre- 
quency, an instability wave develops in the jet shear layer. For clarity, this instability 
wave will be called the fundamental wave or the primary wave. In this work, we will 
use the words “primary” and “fundamental” interchangeably. During the initial stages 
of development, where the amplitude of the disturbance is still small, the stability char- 
acteristics of this wave are described quite accurately by the primary stability analysis 
developed in Chapter 2. However, as this wave propagates downstream, its amplitude 
grows and the mean flow diverges. The mean flow divergence results in the reduction of 
the mean flow vorticity and hence, contributes to a decrease in the local growth rate of 
this wave. Farther downstream, at approximately one or two diameters from the jet noz- 
zle, this wave reaches a finite-amplitude (i.e., saturation) as a result of nonlinear effects. 
The flow field arising from the presence of the fundamental wave may (roughly) be viewed 
as convecting discrete vortex rings. This flow is periodic (approximately) in time and in 
space. 

A principal objective of this work is to determine the sensitivity of this periodic flow 
to various controllable jet parameters such as wavenumber of excitation, azimuthal mode 
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number, Mach number and jet temperature ratio. 

The investigation is carried out via a secondary instability analysis . The instability 
analysis of small disturbances — which for clarity, will be called secondary disturbances 
— superimposed on a base flow consisting of the finite-amplitude fundamental wave and 
a long time averaged mean profile is termed secondary instability analysis. In order to 
make the secondary instability analysis as simple as possible, we invoke the following 
assumptions: 

1. The locally parallel mean profile and the shape of the fundamental mode is minimally 
affected by the nonlinear self- interactions of the fundamental (i.e., shape assumption). 
In physical terms, this fundamental wave can be thought of as a coherent structure 
or a large scale instability wave. 

2. The amplitude A of the fundamental wave, subject to an appropriate normalization, 
remains approximately constant during the evolution of the secondary instability 
disturbance. 

3. The mean flow remains unchanged and the fundamental mode keeps its (initial) shape 
and phase velocity during the evolution of the secondary instability disturbance. 

4. The magnitude of the secondary disturbance is kept sufficiently small to ensure that 
a linear secondary instability theory is valid. 

The justification of the first assumption has been discussed in the primary stability anal- 
ysis. Because of the nonlinear terms in the governing equations, a single linear instability 
mode can never exist, by itself, in the flow. It will interact with itself to produce higher 
harmonics and a mean flow distortion. The effects of the higher harmonics and of the mean 
flow distortion in subsonic round jets were recently investigated in detail by Jarrah (1989). 
His results indicate that the distortion of the mean flow due to wave Reynolds stresses is 
quite small locally and is probably unimportant for the secondary instability. The small 
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change in A during the course of its nonlinear evolution relative to an anticipated strong 
convective growth of a subharmonic disturbance justifies the use of assumption (2). This 
assumption is further supported by experimental studies on planar mixing layers; for ex- 
ample, Ho As Huang (1982) found that when discrete vortices develop on the flow, the 
subharmonic mode grows rapidly while the amplitude A remains approximately constant. 
The recent work of Petersen As Samet (1988) supports assumption (3), especially on the 
insensitivity of the mode shape to nonlinearities, although no detailed experimental results 
exist for heated and compressible round jets in this regard. Furthermore, in subsequent 
analysis, we find that the growth rate of the subharmonic wave is indeed very much larger 
than that of the primary wave; this finding will, in fact, justify the use of assumptions (2) 
and (3). 

In order to have confidence in our results, two different approaches for secondary in- 
stability analysis are used here. The first is based on the method of multiple scales, while 
the second is done by generalizing the normal mode analysis from a parallel base flow to 
a periodic base flow. 


3.1 Method of Multiple Scales 

Kelly (1967) was the first to use the method of multiple scales to describe the- 
oretically the strong growth of a subharmonic disturbance in an incompressible planar 
mixing layer. His results are consistent with the the experimental observations of the 
vortex-pairing phenomenon found by Sato (1956). By using full numerical simulations on 
the instability of an array of Stuart vortices in a mixing layer, Pierrehumbert Widnall 
(1982) confirmed that the two-dimensional subharmonic modes responsible for the coales- 
cence of vortices are indeed the most dominant instability modes. Recently, Monkewitz 
(1988) generalized Kelly*s work and examined vortex pairing in more detail in incompress- 
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ible mixing layers. 

The formulation of the secondary instability problem is based on Kelly’s ideas, although 
a compact notation is used here. The finite- amplitude of the fundamental wave is regarded 
as a moderately small perturbation parameter 8 [say, 8 = 0{A) ^ 0.1]. The physical 
variable, which describes the perturbations, F, is assumed to have a general asymptotic 
expansion of the form 

F = 8Yi + e¥ 2 + 8 2 F n + e6F 12 + 0(e 2 ), 0 < e < 8 « 1 (3.1) 

where Fi is the fundamental (or primary) wave, F 2 is the secondary disturbance (which 
will later be taken as a subharmonic disturbance), Fn denotes the quadratic nonlinear 
interactions of the fundamental wave with itself, and F 12 represents all the quadratic 
interactions between Fi and F 2 . As the Fi mode propagates downstream, it travels into 
a shear flow whose transverse length scale continuously increases, the growth rate of this 
wave continuously decreases and the wave eventually becomes neutral. In shear layers, a 
strong interaction between the fundamental wave and its subharmonic (i.e., one half the 
fundamental frequency) occurs where the fundamental becomes neutral and its amplitude 
is a constant. This neutral state, which arises from the spreading of the shear layer and 
nonlinear effects, is really of no importance in this study. What really matters is that 
there is a region in the jet where the amplitude of the fundamental (or primary) wave is a 
constant and the growth rate of the subharmonic is enhanced from that given by parallel 
flow stability theory; of course, this amplitude is “unknown” to us, although experimental 
data suggests that it can be as large as 10% of the mean flow, and the cross-space structure 
of this wave may be represented by the local linear instability mode of the jet (Petersen 
Samet 1988). In this region, the primary wave interacts with itself at 0(8 2 ) and higher. 
This interaction is assumed to have a negligible effect on the mean flow and the primary 
wave. Consequently, the Fn term in (3.1) can be ignored. Of course, the nonlinear 
interactions of primary wave will eventually lead to the appearance of higher harmonics, 
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but they play no significant role in the generation of the subharmonic at 0(e5) because 
they occur at 0(5 2 e) and higher . 

When the primary wave saturates, the characteristic magnitudes of the primary and 
secondary disturbances are 5 and e respectively. The focus of this work lies at order 
e6 where the primary wave Fi interacts with the subharmonic disturbance F 2 to form 
sum and difference modes. If the axial wavenumbers, frequencies, and azimuthal mode 
numbers satisfy the so-called resonance conditions, which will be discussed momentarily, 
then the growth rate of the disturbance F 2 may be enhanced considerably from that given 
by linear parallel flow analysis. 


3.2 Kinematic Resonance Conditions 


We now illuminate the physics governing resonant interactions between two linear 
waves of the form 

Fi = fi(r)£i + c.c. (3-2) 

F 2 = f 2 (r)i£ 2 + c *c. (3-3) 

where 

E\ = exp * (a\x + exp a\ t 
E 2 = expt(a 2 a: + /? 2 <f>) exp a 2 t 

and c.c. denotes the complex conjugates of the preceding terms. 

After substituting (3.1) - (3.3) into (2.17) - (2.19) and equating coefficients at various 
orders to zero, we obtain the corresponding disturbance equations. The zeroth order 
solutions come from the 0(5) and 0(e) terms. They are the eigenmodes of linear parallel 
flow stability analysis (2.33), and consequently need no further exploration. 
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The equation which governs the quadratic nonlinear interaction between two linear 
waves, obtained at 0(c£), is of the form 

Li 2 (F)i 2 = 22+R+ + 22— R- + c.c. (3 .4) 

where 

E+ = exp[i(ai + a 2 )x + (ft + ft)<£] exp(<Ji + <7 2 )t (3.5) 

E- = exp[i(aq - a 2 )x + (ft - ft)$exp(<7i + <r 2 )t , (3.6) 

R+ and R_ are given in Appendix C, L 12 is the linearized operator (given in Appendix 

A) which governs the evolution of small disturbances in a parallel base flow, the tilde 

denotes the complex conjugate of a quantity, and c.c. is the complex conjugate of the 
terms preceding c.c.. 

Since L 12 is a linear operator, the solution F 12 consists of the superposition of the 
homogeneous solution and the particular solution due to the forcing function, namely the 
right hand side of (3.4). The homogeneous problem satisfies 

1*12 (F l2)hom — 0 (3.7) 

and the particular problem is determined from 

2^12 (F \ 2 )pat — 22+R+ + 22_R_ + c.c. (3.8) 

where the subscripts pat and hom respectively denote the particular and the homogeneous 
solutions of the inhomogeneous problem (3.4). Therefore, the physical variables can be 
written, up to order eS y as 

F = SFi + el? 2 + e6[{Fn)hom + (Fn) P at) (3.9) 

Since the homogeneous solution (Fn^om has an initial magnitude of eS and is described 
by linear eigenmode analysis, it will always continue to be smaller than F 2 by an order 5, 
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and hence, can be neglected. The interactions between the two waves of the form given 
by (3.2) and (3.3) will produce the following four sets: 

(a\ + + 02 ><rx + ^ 2 ), {-ot\ + <*2, —fix + ^ 2,^1 + 02), 

(Ofl - OL 2 ,fix - fii.vx + ^ 2 ), (-oti - <X2,-fix - fit, Zx + h) (3.10) 

As long a is these four sets do not satisfy the dispersion relation (a,/?,cr), a solution to 
(3.4) may be found, in principle, in a straightforward manner by the method of variation 
of parameters. In such cases, (3.4) represents am inhomogeneous equation for which a 
particular solution may be sought in the form 

(Fu)pa* = £+F2i(r) + E-V^r) + c.c. (3.11) 

Equation (3.11) shows that the growth or decay of (Fi2)pat is the same as that given by 
linear parallel flow analysis. Further, in order to satisfy the resonance conditions, which 
will be discussed momentarily, the disturbance F 2 is required to be the most unstable wave 
given by parallel flow theory. This implies that, for the case where resonance conditions 
are not met, the (Fi2) pa t will remain smaller than F 2 by 0(5), and it can therefore be 
similarly neglected. 

On the other hand, when one of the four sets of a, fi> and a satisfies the dispersion 
relation, the operator L 12 , with the appropriate boundary conditions, admits nontrivial 
solutions which are, in fact, the linear instability modes of the base flow. This leads 
to the appearance of secular terms in (3.4) which results in the occurrence of resonance 
interactions. Accordingly, a solution is possible only if the right hand side of (3.4), which 
contains the resonance terms, satisfies the so-called solvability condition (or Fredholm 
alternative). In general, this constraint cannot be satisfied using a regular perturbation 
expansion. Before we proceed to consider how this constraint can be determined from a 
singular perturbation method, let us now explore the resonance interactions a step deeper. 
To do this, we consider the linearized in viscid primary disturbance equation (2,48) and its 
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conjugate; they are, respectively, 


s*<;- 


K 

«o — /ot 
Uo + iff /a 


-&)%■+ [« J (»o - CfeoMj a ! ]p = 0 

Po dr r* 

~^)f+ [a*(«o - Cf»M/ a *|p = 0 

Po dr 4 r* 


(3.12) 

(3.13) 


where C is the complex phase velocity, <r — — taC, the primes denote differentiation with 


respect to r, and the tilde denotes the complex conjugate of a quantity. The complex 
phase velocity of an unstable mode which satisfies (3.12) and (3.13) can be expressed as 


C(a,£) = C*(|a|, |/J|) + sjn(a)»'C/(|a|, |/J|), (3.14) 


and the corresponding dispersion relation is 

ff(aj) = MC/(|a|, \P\) - iaC R {\a\, |/?|) (3.15) 

where C = Cr + tCj, and s^n(a) means the sign of a. In view of (3.10), the complex 
equation (3.15) can be separated into two real equations — a growth rate equation and a 
frequency equation — which govern the resonance conditions between two unstable waves . 
These equations are, respectively, 

|ai|C/(M, \Pi\) + |«2|C f /(|«2|, \P*\) = | <*2 + a i|£/(l a 2 + «i|,|A + &I) (3.16) 

aiCjz(|oi|j |£i|) + oiOjtdajUAI) = (a2 + <*i)Cj?(|a2 + <*iUA + £ 2 1) (3.17) 

where the a’s and s can be positive or negative, and are so far arbitrary. It may be 
noted that when /9i = 0, the resonance conditions are valid independent of the direction 
to which a wave is propagating; this is because the inviscid equations (3.14) and (3.15) do 
not distinguish /? from — /?. 

In order to satisfy the resonance conditions, the wavenumbers and mode numbers of 
two unstable waves are restricted to the following choices: 


1. The fundamental mode is axisymmetric and has an axial wavenumber close to a 
neutral wavenumber, [i.e., /?i = 0, C/(|ai|,0) ^ 0]. 
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2. The secondary disturbance, whether axisymmetric or helical, is an almost maximally 
amplified subharmonic wave, [a* — ai/2]. 

It may be remarked that although choices (1) and (2) are necessary, they do not guarantee 
resonance because when the waves are dispersive, condition (2) places a restriction on the 
phase velocity, Cr. 

When the resonance conditions are (nearly) met, a subharmonic wave interacts effectively 
with the fundamental to produce an instability wave whose complex growth rate, < 7 , is 
(nearly) the same as that of the subharmonic. It is effective because resonant interac- 
tions produce secular terms which lead to a particular solution of (3.4) that grows more 
rapidly than the corresponding homogeneous solution by at least a factor t. Of course, 
the wavenumber and the mode number of the wave produced is always the sum of the 
wavenumbers and mode numbers of the two interacting waves. Consequently, an axisym- 
metric subharmonic which interacts with its fundamental will reproduce itself. On the 
other hand, a helical subharmonic mode which interacts with the fundamental will not 
reproduce itself but will produce a mode that propagates at an equal but opposite angle 
to the z-direction. For example, an axisymmetric mode with the wavenumber 2a and' 
the mode number 0 — denoted by (2a, 0) — interacts with a helical mode (—a, fl) to 
excite a mode (a, /?). Similarly, the mode (2a, 0) interacts with mode (-a, —/3) results in 
the mode (a, — /?). Therefore, an interaction between mode (2a, 0) with modes (-a,±/?) 
results in modes (a,±/3). Because a helical mode which interacts with the fundamental 
mode does not reproduce itself, it is necessary to consider the secondary disturbance to 
be comprised of a pair of helical waves spinning at an equal but opposite direction to the 
x-axis. Consequently, if a non-axisymmetric subharmonic is involved, we need to assume 
that the secondary disturbance is of the form 


F 2 = ^2r + + C.c. 


(3.18) 
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where 

F 2r = f2r(0 exp t(a 2 X + /? 2 <£)expa 2r * + c.c. (3.19) 

F 21 = f 2 j(r) expi(— a 2 x + /? 2 <£)exp<7 2 j£ + c.c. , (3.20) 

where the subscripts r and / denote a helical wave propagating in the direction of positive 
and negative <f> y respectively; for clarity, we call these waves right and left spinning (or 

propagating) waves, respectively. Since the terms involving the product of e are trun- 

cated in the expansion (3.1), the left and right spinning waves of F 2 do not interact with 
each other. Each will, however, independently interact with the primary wave Fi- For 
example, F 21 interacts with Fi, producing F 2r , and F 2r interacts with Fi, producing F 2 j. 
Consequently, F 2 is again reproduced. Furthermore, since both the mean flow and the 
axisymmetric Fi are independent of the azimuthal angle, we can assume that the growth 
rate of the right and left spinning waves, as a result of resonant interactions between F 2 
(given by 3.18) and an axisymmetric primary wave, are modified exactly the same. In 
other words, there is no physical preference for the growth rates between the left and the 
right spinning waves. In fact, a second approach — which will be discussed in section 3.5 
— shows that the growth rates of F 2r and F 2 j in a frame of reference that moves with the 
phase velocity of Fi are indeed the same. 

When the resonance conditions are met, the subharmonic disturbance interacts strongly 
with the fundamental to reproduce itself and therefore to increase its growth rate. The 
modified subharmonic growth rate, as a result of resonant wave interactions, will be ac- 
commodated by a set of slow scales. 

We now proceed to obtain the governing equation for the modified growth rate using 
the method of multiple scales. The method of multiple scales is discussed in several books 
[e.g., Nayfeh 1973]. Accordingly, we can write 

F = <5Fi + eF 2 + e8Ti2 + 0 ($ 2 ), 0 < c « 8 <£ 1 (3.21) 


where 
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Fi = Fi(x ,t,<j>,r) 

(3.22) 

F} = F 2 (x,n,t,ti,<f>,<f>i,r) 

(3.23) 

r 12 = F ii(x,xi,t,ti,<f>,<l>i,r ) 

(3.24) 


and = 8t,<f>\ = 8<f> y x i = 8x are the slow scales of time, azimuthal angle, and axial 
direction, respectively. It may be remarked that although this set of slow scales eliminates 
the nonuniformity caused by the secular terms in equation (3.4), it does not determine a 
uniform zeroth order expansion. This is because the amplitudes of the mean flow and the 
fundamental wave are assumed to remain unchanged with time, while the subharmonic 
wave, chosen from the resonance conditions, is strongly amplified. The series expansion 
(3.21) becomes disordered whenever the total magnitude of the growing subharmonic 
exceeds the preceding term; uniform validity only holds, at most, on a time scale 8~ l . 
This set of slow scales is compatible with (3.21) where the 8 2 term which arises from the 
third order interactions (i.e., 8 2 e) is truncated. It is used to exam i n e the influence of the 
fundamental wave on its subharmonic to linear order in terms of the amplitude 8. In the 
following equation (3.27), we shall assume that the shape of the primary subharmonic is 
unaltered by the appearance of the fundamental. This assumption is more fully justified 
by the second approach, to be discussed in the next section, which allows the shape of the 
subharmonic in the presence of its fundamental to be determined. 

Since the physical variable F is a function of fast and slow scales, the partial derivatives 
in t, x, and <f> become 

dt dt "I" dt\ 

dx dx dx\ 

A 8 — 

d<j> d<t> + d d<f>!' 


(3.25) 
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As in the parallel flow analysis, we assume the normal mode concept and write 

Fi = fi(r)£i + c.c. (3-26) 

F 2 = fl(*i,*i,^i)f2(r)^ + c.c. (3 27) 

where 

Ei = exp»(aix + exp ait 

E-i = exp «(a 2 x + fo(f>) exp <r 2 t 

and c.c. denotes complex conjugate of the preceding terms. It should be noted that in 
(3.27), we have omitted the left propagating helical component with the understanding 
that the left propagating helical wave which interacts with an axisymmetric primary wave 
produces and, therefore, enhances the growth of a right propagating helical wave; similarly, 
the left propagating helical wave can be reproduced through the interaction between a right 
propagating helical wave and an axisymmetric wave. 

The amplitude function B is slowly varying in space and time, and will be determined 
at 0(e£) by means of a solvability condition; the mode shape f 2 is that given by the parallel 
flow theory. 

It is now perfectly clear from equation (3.27) how the growth rate of the subharmonic is 
modified. Because of the departure of the growth rate of F 2 from the parallel flow analysis, 
a set of slow scales is needed to reflect this change. For example, the slow scale Sx implies 
a long length scale of 0(R*/6), in addition to the usual length scale R* . This means that 
as long as the evolution of the subharmonic is on a length scale that is much smaller than 
0(R* /S), its local growth is given by parallel flow theory. However, if the growth of the 
disturbance is on a distance that is comparable to 0(R* /S), the change in the amplitude 
function B will then be of the same order as that provided by the parallel flow analysis. 
In temporal theory, it is clear that this change, which may be interpreted as an additional 
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(positive or negative) growth rate for the disturbance, is physically meaningful only over 
a time scale that is of 0 ( 1 ). 


3.3 Amplitude Equation 

We now proceed to derive the amplitude equation of the subharmonic disturbance 
when the resonance conditions are met. The amplitude equation will be derived using 
two independent methods: the first is based on an “improved” inviscid analysis and the 
second is based on a viscous theory. These derived amplitude equations will be called 
an inviscid amplitude equation and a viscous amplitude equation, respectively. Here, the 
word “improved” signifies that the primary instability modes are actually obtained from 
viscous considerations. The inviscid amplitude equation is derived from a single ODE, 
while the viscous counterpart is based on a system of ODE’s. Since the wavenumber of 
a subharmonic wave is near a maximally amplified wavenumber, the growth rate of the 
subharmonic obtained from an inviscid analysis will provide a good approximation to that 
obtained from a fully viscous analysis. We now proceed to derive the inviscid amplitude 
equation, while the viscous amplitude equation is given in Appendix D. 

For the particular case where the product of Fi and F2 reinforces F2, the disturbance, 
F12, can assume the form: 

F 12 = + c.c. ~\ ( 3 . 28 ) 

where 

^12 = (« 12 , VtfyWUjPUyPu) 7 , 

E 2 = exp t (0:2 x + j3 2 <j> ) exp a 2 t 

and the dots stand for the non-resonating terms that are unimportant. 
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We note that the linearized primary disturbance equations for an inviscid compressible 
jet at 0(e8) can be manipulated to a single pressure perturbation pi 2 (r) of the form 


Mpu) = ^r + 


1 _ 2 Up _ Pq' dpi 2 

r uq — C 2 po J dr 


fil 


a 2 («o - C 2 ) 2 PoMy - -i - «2 


p 12 = rhs 

(3.29) 


where 


as f as , as , _ & 

r/is — hi h hs^T” + hi2# 

ati dzi <?<pi 

2tiQ Pq d 


(3.30) 


— — + — ](poV2) — i^2P0^2/ r + *^ 3(^0 “ C2)PoMj P2 


. r X 

kl = ~ ta ^2 " I" - , 0 ' * 

h 2 = -m 2 p 0 u ou 2 - [i - Jfs c ~ ^ + ^J( UoPot '^ " 

+ia 2 (uo - C 2 )poAd/(«oP 2 + 7Po« 2) 
ha = ia(uo — C^poUvi/r + M? r P 2 ) 


hi 2 = ia 2 * tt - + [i - -^2- - £° + i-ln- + ifiiFv>-/ r ~ i* 2 («o - C 2 ) P0 MjF p - 
1 r tio — C2 Po <* r 

In the above equations, the tilde denotes the complex conjugate of a quantity, the prime 
denotes a differentiation with respect to r, C 2 is the complex phase velocity of the subhar- 
monic mode, while F u -, F v -, K- and F p~ are 8 iven “ Appendix E. Note that although 
(3.30) is derived based on an inviscid analysis, its coefficients, h's, are evaluated using 
viscous parallel flow instability theory. Further, h 12 represents the particular set of the 
interaction terms between F x and F 2 that give rise to a subharmonic disturbance, while 
hi, h 2 , and h 3 come from the dependence of F 2 on the slow scales. 

We multiply (3.29) by r/{po(«o - C 2 ) 2 } to make it self-adjoint. The equation in the 

self-adjoint form becomes 

C(pu) = — (t»o — Ci) 2 L(pu) = rhs ( — («o — C 2 ) ) 


(3.31) 
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In order for (3.31) to possess a nontrivial solution, its right hand side must be orthogo- 
nal to all i’s such that £z = 0 where £ is the adjoint operator and z is the adjoint solution. 
The orthogonality requirement, or solvability condition, is necessary and sufficient for the 
existence of solutions at order 6e. We note that since the operator £ is self-adjoint [i.e, L 
— £] and the boundary conditions on z are the same as those of the homogeneous solution 
[i.e., C(pit) = 0], the adjoint solutions are identical with the homogeneous solution. 

The orthogonality requirement becomes 

dB dB dB ~ „ 

+ g 2 — + g 3 — + ?4 S = °. (3.32) 

where 


f°° r ,1 2t4 p' Q . 

<Il= Pi 7 77T7 {*«2POW2 + ( -p; - —)PO v i 

Jo Po(«o — Ci) 2 r uo — C po 

+iPipow 2 /r - i'a 2 (uo - C 2 )poMjp 2 } - ( c 7 ) 2 ^ p ° V2dr 




> r f . * ,1 2tiQ 

P 2 7 — > 2 + ( 77 

~ Cl) r u o "" C2 Po 




+i0iK-/r ~ «'« 2(1*0 - C 2 )p Q MfF p -} - { —^-^ )'F v -dr, 

where the (•)' denotes a differentiation of its argument with respect to r, and q* con- 
tains the quadratic resonant wave interactions. Here q 2 and <73 are included to allow the 
resonance to be manifested spatially. Their explicit expressions are not given here because 
this thesis focuses only on temporally growing waves. The integrands in the q’s are eval- 
uated from linear, viscous, parallel flow instability analysis. It may be noted that if we 
include the Fn term in the series expansion (3.21), an amplitude equation for the modi- 
fied growth rate of the fundamental can similarly be obtained at order <5 2 (see Monkewitz 
1988). This modified growth rate, which will be coupled into equation (3.32), provides a 
relationship for the mutual feedback between the subharmonic and the fundamental. In 


accordance with our previous discussion, we shall not obtain the amplitude equation for 
the primary wave because this amplitude is assumed to be a constant. The constant ampli- 
tude A assumes a catalytic role in the enhancement of the growth rate of an infinitesimally 


small subharmonic secondary disturbance. Therefore, we called our instability analysis 
a linear “secondary” instability analysis instead of a “weakly” nonlinear analysis for two 
interacting waves. 

The equation (3.32) together with its complex conjugate yields 


, l2 d 2 . ,2 d 2 . l2 d 2 

^ dt\ + ^ dx\ + ^ d<f>\ + 

d 2 d 2 d 2 

Real[2$ig 2 sm + ^7 + 2 g2g3 a ‘ = 


' dt\dx\ 


dtid<f> 1 




(3.33) 


where the tilde stands for the complex conjugate of a quantity. The above equation is 
valid for temporal or spatial analysis because the amplitude B is a function of spatial and 
temporal variables. Since the coefficients in (3.33) are independent of the slow scales, we 
can write 


B ~ exp(Aixi + Ajii + As^i) 


(3.34) 


and substituting the above into (3.33) yields 
|gi| 2 A* + ]? 2 | 2 A 2 + |gs| 2 Ai + (2§ig2)i?AiA2 + (2gig3)jeA2As + (2g2?3)flAiA3 = |g4 | 2 (3.35) 


The A’s are in general complex and ( • )r denotes the real part of a complex number. 
Each real part of Ai, A 2 , A 3 modifies the streamwise, temporal and azimuthal growth rate, 
respectively. Similarly, each imaginary part of A’s modifies its corresponding frequency (or 
wavenumber). Running parallel to Kelly’s investigation of the mixing layer, our investiga- 
tion is based on a disturbance of a fixed spatial wavenumber growing in time. Therefore, 
Ai and A 3 are zero, while A 2 is real. (3.35) can therefore be simplified to 
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The secondary subharmonic disturbance can be expressed in terms of the original time 
variable as 


Fj = fa expt*(a 2 £ + <f>) exp{<72 + A 2 $)£ (3.37) 

The dependence of this growth rate as a function of various jet parameters is discussed in 
Chapter 6 . Clearly, the growth rate of the subharmonic wave is changed by A 2 6 due to 
parametric resonance. 


3.4 Normal Mode Analysis 

Our previous approach models the secondary instabilities of a round jet, which 
results in vortex pairing, by a perturbation analysis of two instability waves, Fj and F 3 , 
of different wave vectors undergoing interactions. When their product reproduces with 
F 2 , the growth rate of Fj is modified. While this model provides the physics of resonant 
interactions, it has at least three shortcomings. 

1 . The finite amplitude of the fundamental wave is regarded as a small perturbation 
parameter. Therefore the perturbation theory, can provide an accurate solution only 
for those values of 8 for which the perturbation series is an accurate representation 
of the exact solution. As usual, the region of validity of an asymptotic analysis is 
not known a priori, and remains to be determined from a more complete analysis. 

2. There is no mechanism for the determination of the threshold amplitude of the 
fundamental where the phase- locking with the subharmonic occurs. It is now known 
that this phase-locking is a mechanism for the onset of secondary instability. 

3. As previously mentioned, this perturbation approach fails to determine the correction 
to the shape of a secondary subharmonic mode; it only provides information on its 



66 


growth rate. 

These shortcomings can be circumvented by an alternative approach developed by Her- 
bert (1983). He pioneered the development of the normal mode analysis of the secondary 
instability arising from the presence of TS waves in plane channel flows. His theory pre- 
dicts with reasonable accuracy the growth rate of a 3-d secondary disturbance up to the 
stages prior to laminar breakdown and transition. In his subsequent analysis of the (lin- 
ear) secondary instability of flat plate boundary layers, he has shown that this model is 
capable of explaining all the essential flow physics of staggered lambda vortex formation 
(1984) and peak-valley splitting (1985) in Blasius flows. 

Motivated by Herbert’s success, we will adapt his approach to subsonic heated round 
jets. In this approach, a sufficiently small disturbance F 2 is superimposed on a periodic 
base flow consisting of the steady unidirectional mean flow Fo and a finite amplitude 
axisymmetric primary fundamental wave Fi. We assume that the mean flow as well as 
the shape of Fi is minimally affected by the nonlinear self-interactions of Fi. The choice of 
an axisymmetric wave is based on the resonance conditions established previously. Using 
a frame of reference moving downstream with the phase velocity of Fi, it is appropriate 
to decompose the perturbation vector field F as 

F (x,r,<f>,t) = AFi(x,r) + eF 2 (x,r,<f>,t) (3.38) 

where 

Fi = f\(r)e iax + c.c. , (3.39) 

A y as discussed earlier, is assumed to be a constant, and x now refers to the streamwise 
coordinate of the moving frame of reference. In (3.38), A is usually small but finite (say, 
A ^ 0.1), not a perturbation parameter, while e is assumed to be sufficiently small to allow 
linearization about a steady and spatially periodic base flow. Of course, the range of e for 
which linear theory is valid remains to be determined by comparing it with future work, 
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say, direct numerical simulations. The finite value of A will remove the usual accuracy 
issue associated with an asymptotic analysis. Further, since the amplitude A has a physical 
meaning provided the eigenfunction fi is suitably normalized, we choose, for convenience, 

max |ui(r)| = l/2 (3.40) 

0<r<oo 

so that A is the maximum magnitude of the streamwise velocity perturbation. This nor- 
malization, which is also used in our multiple scales approach, is not unique and is rather 
subjective because fi is an eigenfunction. Herbert normalized A so that it is the maxi- 
mum rms value of the streamwise velocity fluctuation. The magnitudes of eigenfunctions 
of Cohen (1986) are normalized such that the experimental and theoretical profiles have 
identical areas. At this point, it is important to note that this secondary instability anal- 
ysis is still linear, as is the classical instability analysis. In the latter, the flow is linearized 
about a locally parallel flow . The coefficients in the resultant equations depend solely on r 
since the (x, <ft) and t dependences of the (primary) disturbance are trivially extracted by 
Fourier-Laplace transforms. This leads to the mode shape fi(r). In the present (secondary 
instability) analysis, linearization is made about a flow which is periodic in the streamwise 
coordinate x, This streamwise periodic flow becomes independent of time in a reference 
frame that moves with the phase velocity of the fundamental wave, and is independent of <f> 
since Fi is axisymmetric. Therefore, the <f> and t dependence of the secondary disturbance 
can again be extracted by Fourier - Laplace transforms, but the x dependence requires a 
more careful consideration and additional assumptions, to be discussed shortly. We can 
now apply the normal mode analysis in <f> and t and write 

F 2 — ^(x, r)e l/9a ^e* 3 * + c.c. , (3.41) 

where is the complex growth rate, and /? 2 is the azimuthal mode number. After substi- 
tuting (3.38), (3.39) and (3,41) into the governing equations (2.17) - (2.19), and neglecting 
quadratic and higher order terms in e, we obtain a system of linear partial differential equa- 
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tions which govern the instability of the linearized secondary disturbance, F 2 . Apart from 
a complicated r dependence and numerous other terms, these equations are essentially of 
the classical Mathieu type because the coefficients are (streamwise) periodic. It is known 
that, according to Floquet theory of linear periodic differential equations, a Mathieu equa- 
tion whose coefficients are x-periodic has two important classes of Floquet solutions: a 
7 r-periodic and a 27r-periodic function. With guidance from Floquet theory (Ince, 1956), 
we can therefore express the secondary disturbance in the Floquet form as 

h (x,r)->e***i 2 {x,r) (3.42) 

where a x is the characteristic exponent and is periodic in x. The solutions with period Aj 
— where Ai is the wavelength of the primary wave — are called the primary (fundamental) 
resonance, while the solutions with period 2Ai are known as the principal parametric 
(subharmonic) resonance. The instability modes which arise from these resonances can, 
in general, be expressed in terms of Fourier series as 

£ f2 Ar)^ iaix/2 (3-43) 

y=-oo 

The fundamental modes are the harmonic series whose period is the same as the streamwise 
wavelength of the primary wave and are given by j = even. The subharmonic modes have 
odd j, and therefore, constitute a subharmonic series whose basic wavelength is twice 
that of the primary wave. The fundamental and subharmonic modes are completely 
uncoupled since the wavelength of the basic flow is different from the basic wavelength of 
the subharmonic series. In this work, we shall investigate only the parametric resonance 
case. 

We now digress briefly to present some known results in incompressible boundary layer 
flows. The fundamental modes produce the peak-valley splitting of the lambda-shaped 
vortical structures, while the subharmonic modes cause the staggered formation of lambda 
vortices. By approximating the Fourier series (3.43) to consist of merely two oblique waves 
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propagating at equal and opposite angles to the z direction, Herbert (1984) has shown that 
these two oblique subharmonic modes, which provide a reasonably accurate approximation 
to the streamwise structure of the subharmonic disturbance, are sufficient to give rise to 
the formation of staggered lambda vortices. 

We return to equation (3.43). Its substitution into the linearized disturbance equa- 
tions yields an infinite set of ODE’s for the complex growth rate &2 and the mode shape 
Uj- Numerically, the streamwise structure of the secondary disturbance is expanded in 
terms of a finite Fourier modes whose stability is governed by a system of N a ODE’s — 
where N a is 5 times the number of Fourier modes used. Numerically, these (stability) 
equations are approximated by a complex algebraic eigenvalue problem whose computa- 
tion requires O (N*) storage and O (N*) operations (see Chapter 4) — where N i equals 
N a times the number of grid points in the computation domain. In order to be practi- 
cal in the numerical evaluation of the secondary stability problem and in order to make 
a direct comparison with the previous perturbation approach, simplification of equation 
(3.43) is deemed desirable. The justification of this simplification will now be discussed. 

In controlled jet experiments, a primary instability wave, which is often generated by 
an acoustic source or through vibrating ribbons, produces a vortex structure that is re- 
sponsible, at leasj^ partially, for the mixing and spreading of jet shear layers. In order to 
control the global spreading rate via the manipulation of this vortex structure, a discrete 
subharmonic disturbance can be excited simultaneously with the primary (or fundamen- 
tal) frequency. In accordance with our discussion in section 3.3, when the subharmonic 
disturbance is axisymmetric, it is assumed to be a subharmonic wave; when the subhar- 
monic is helical, it is comprised of a pair of helical waves whose propagation angles are 
equal but opposite with respect to the mean flow direction. As dictated by the reso- 
nance conditions, these excited discrete modes are sufficient to promote a strong growth 
of a subharmonic disturbance. In order to compare the two approaches, we assume that 
the (secondary) subharmonic growth rate is minimally affected by the unexcited higher 
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terms in the subharmonic series. In other words, we assume that the (two) Fourier modes 
(ax/2 and — ai/2) are reasonably sufficient to represent the the streamwise structure of 
the subharmonic disturbance. 

With the preceding assumptions, we can approximate the subharmonic disturbance by 
F 2 = V J V^(f 2(1 e’ a */ J + f 2 ,-ie- iai/2 ) + c.c . , (3.44) 


where 


h t i = (U2 f l, t/2,1, U>2,1 ,P2 t l>P2 i l), 


and, 

*2,-1 = (^2,-Ij V2 ( -l, W2,-UP2,-liP2 t -l) 

As in primary stability analysis, we employ temporal theory, which implies that fr x = 0 
and &2 ^ 0 is the complex eigenvalue to be solved. Thus, in a temporal theory, we can 
write 

F 2 = e ^V' ? 3 *(f 2 l le , ' a, / 2 + f 2 | -ie _,ai/2 ) + C.c. (3.45) 

Note that in the moving frame of reference, (3.45) shows that the complex growth rates 
of a pair of helical modes, whose propagation angles (with respect to the z-axis) are 
equal and opposite, are exactly the same. 

To recall, a main result of the multiple scales analysis is that the secondary disturbance 
can interact with its fundamental to produce a wave whose wavenumber and frequency are 
the same as the disturbance itself and can therefore reinforce it. Since such a destabiliza- 
tion of the subharmonic requires the matching of its phase velocity with the fundamental, 
our calculations focus on the particular case where the subharmonic becomes phase-locked 
with the fundamental and <r 2 becomes purely real. The modal representation of F 2 , given 
by (3.45), will transform the linearized secondary disturbance equations from a system 
of partial differential equations into a system of ODE’s. These ODE’s are obtained after 
substituting (3.38), (3.39) and (3.45) into (2.17) - (2.19), and by neglecting quadratic and 
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higher order terms in e. The resulting equations become 

(A 0 + AAi + A a A 2 + A s A s )g = * 2 fig (3.46) 

where 

g = (u2,l> t>2,— 1 j ^2,1 , **>2,-1, P2 t l>P2,-l> P2,U P2,-l) T ( 3 -47) 

Recall that when the secondary disturbance is axisymmetric, we need only to consider 
one subharmonic wave. In this case, f 2 ,-i = f 2 ,i where tilde is the complex conjugate of a 
quantity. In (3.46) the precise elements of Aq, Ai and B can be extracted from Appendix 
F, the A 2 term arises from the cubic interaction, A 2 e , while the A 5 term originates from 
the quartic interaction, A s e. The quartic term stems from the V 2 (p/p) term of the energy 
equation; of course, this term would have disappeared if the temperature had been used 
as a dependent variable. Note that for incompressible and unheated flows, A 2 and A 3 do 
not exist because the nonlinearity in the governing equations is quadratic. Further, since 
our representation of the flow arising from the finite-amplitude saturation of the primary 
wave is approximated only to order A (see Chapter 6), we ignored, for consistency, the A 2 
and A 3 terms. Hence, the secondary disturbance equations become 

(A 0 + AAi)g = <r 2 6g (3.48) 

We should recognize that when A ^ 0 in (3,48), the parallel jet mean flow is slightly 
distorted by the presence of the fundamental mode. This leads to a modified growth rate 
of the subharmonic from that obtained through a parallel flow analysis. 

At this point, it is crucial to realize that the difference between our two different 
approaches lies in the treatment of the terms representing the interactions between two 
waves. When using the normal mode approach, the interaction terms appear explicitly as 
the coefficients in the disturbance equations (3.48). In contrast, these terms, when using 
the method of multiple scales, appear as inhomogeneous terms in the equation (3.4). 
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The boundary conditions at r = 0 and r = oo for the secondary disturbances are 
identical with those of the primary stability analysis. The system of ODE’s (3.48) together 
with the appropriate boundary conditions constitute an eigenvalue problem that yields, 
for a given mean flow and primary eigenmodes, a dispersion relation of the form 

o-i - £ 2 ( 0 : 2 , £ 2 ) (3.49) 

The evaluation of the eigenvalues and their corresponding eigenfunctions is described in 
Chapter 4, and the secondary instability results are discussed in Chapter 6. 
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Chapter 4 

NUMERICAL METHODS 


This section deals with the numerical formulation of the eigenvalue problem which 
arises from the analysis of the stability of subsonic round jets. The primary stability 
problem can be cast in a system of 5 coupled ODE (2.33), while the normal mode analysis 
of the secondary (subharmonic) instability problem is approximated by a system of 10 
coupled ODE (3,48). The numerical solution of these systems of equations is obtained by 
a spectral collocation method with Chebyshev polynomials, and the primary instability 
problem is checked by a finite difference method. For both methods, a function, say f[z ), 
is approximated by a discrete set of values of /(z) on a set of nodal points z,. The discrete 
pointwise approximation of the function /(z) will eventually convert a linear differential 
equation into a system of algebraic equations which can be solved by several techniques. 
We begin by outlining an application of spectral method to the instabilities of shear layer 
generated from a jet emanating from an axisymmetric nozzle. 


4.1 Spectral Method 

The spectral method has now been successfully implemented in many types of 
hydrodynamic stability problems such as boundary layers, wakes and mixing layers. A 
detailed and extensive discussion of the theory and application of spectral methods, partic- 
ularly with regard to fluid dynamic problems, can be found in the monograph of Canuto, 
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Hussaini, Quarteroni, and Zang (1987). The starting point of a spectral method is to 
expand a function u = u(z) in an infinite complete set of orthogonal functions of the form 

u(z) = f^a k T k {z) (4.1) 

k=0 

where the T*(z) are the known basis functions and the coefficients a* ( k = 0, 1,* • •) are 
called the spectrum of u(z). Numerically, the function u(z ) is expanded in terms of a 

truncated set of the basis functions T*. For example, 

M 

u ps t ij/(z) = ]T^a*T*(z) ( 42 ) 

Jt=0 

where tijy(z) is an approximation to u(z) for each M > 0. For a smooth function with 
periodic boundary conditions, the natural choice of the basis functions is sines and cosines. 
This leads to the usual Fourier series and its finite version. On the other hand, for non 
periodic boundary conditions, the most widely used basis functions are the Chebyshev 
polynomials. In any expansion, the task is either to find the set of values of u at the nodal 
grids through an appropriate interpolation or the spectral coefficients a* (k = 0, 1, • ■ • M). 
Intuitively, when the number of grid points is sufficiently large, they are equivalent to each 
other. 

The Chebyshev polynomial of degree k , T*(z), is defined by 

Tk{z) = cos(A:cos“ 1 z) — 1 < z < 1. (4.3) 

Therefore, Tq = 1 , T x = z, and the higher polynomials can be generated by the recursion 
relation, 

2 zT k = T k+1 {z) + Tjt-^z), k > 1 (4.4) 

Details of various identities, orthogonal properties, and numerical quadratures can be 
found in the book by Gottlieb and Orszag (1977). The Chebyshev polynomials are even 
for even k and odd for odd k, as illustrated in Figures (4.1a) and (4.1b). They are the 
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solutions of the Sturm-Liouville differential equation 

Vl - z 2 ^ (VI - z 2 ) ^ + k 2 T k = 0 (4.5) 


which are bounded at z = — 1 and z = L 

It is known that if a smooth function is expressed in terms of the eigenfunctions of 
a suitable Sturm-Liouville problem, “spectral accuracy” is obtained. This means that 
the spectrum a*, according to this basis function, decays more rapidly than any power 
of k. To be more specific, the Fourier cosine series for an even function f(x) is used for 
illustration. The appropriate expansion is 

oo 

f(x ) = ^ a* cos(jfcx) 0 < x < x (4.6) 

ifc=0 


where 

a k = — f f(x)cos(kx)dz, k = 0,1,2,'" (4.7) 

x- Jo 

If a periodic function /(x) has At — 1 continuous derivatives, and the At** derivative is 
integrable, then, after integrating (4.7) by parts (At — 1) times, we find that 


a * = Jo f {M) W COS (W 2 - kx ) d * 


(4.8) 


where \ — 1 when At is even and 0 otherwise, and the /W is the A<-th derivative of /. It 
is worth noting that the integral in equation (4.8) is bounded by some constant which is 
independent of k but may depend on At. Therefore, the truncation error, say e, using the 
M + 1 modes to represent an infinitely smooth function by (4.2) can be estimated to be 

const ant ( At) 


e < 


~w 


for any At 


(4.9) 


This is because each of the coefficients a k [for k > M] satisfies (4.8) and the neglected 
series expansion for /, 

oo 

]P a k cos (kx), 

k=JJ 
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will be of the order given by the right hand side of (4.9). [More precise error estimates for 
the eigenfunction expansion in terms of Chebyshev polynomials is discussed in Canuto et 
al] . It may be remarked that a finite difference method with M grid points usually has an 
error estimate of M~ p where p is the order of the accuracy of a particular scheme used; 
for example, a second order method has p — 2. Thus, if a function is infinitely smooth, 
then, according to (4.8), the spectrum a* decays faster than any finite power of k y or, 
simply, it decays “spectrally”. An important consequence of this spectral accuracy is that 
accurate numerical results can be obtained by using relatively “few” grid points, and it is 
precisely this accuracy that makes spectral methods an emerging powerful computational 
tool. 

Further, after making the substitution z = cos(x) into (4.5), we find that 

^± + k*T k =0 (4.10) 

so that the Chebyshev polynomials are the usual trigonometric functions of the argument 
kx = A: cos* 1 z. From (4.3) and (4.6), we obtain 

/(*) = £>**(*), (4.11) 

where x = cos _1 (z). Therefore, according to equations (4.11) and (4.1), the Chebyshev 
series has the same convergence properties as the Fourier cosine series, since the former is 
really the latter in disguise. 

If we expand an arbitrary function u in (M + 1) Chebyshev polynomials, (4.2) will 
contain (M + 1) unknown coefficients, ao, ai, In order to evaluate these using 

collocation methods, the approximate expansion of u(z) is satisfied exactly at a set of 
(JV -f 1) points, say zq, z\ y .., zy. These points are called collocation points; they are 
arbitrarily subjected to a rather qualitative constraint. In order to see how to choose 
these collocation points, we go back to (4.8) and note that the required coefficients a* (A; = 
0, 1, M) are given by a trigonometric integral. This integral can be evaluated extremely 
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accurately by a trapezoidal formula because, in fact, this scheme is equivalent to an M- 
point Gaussian quadrature rule. Therefore, an optimum way is to divide the interval 


0 < x < x into equal increments: 


Xj = T j = 


(4.12) 


or in terms of the original variable, z, 


Zj - cos j = 0,1, ...,M 


(4.13) 


The set {zj} is called the Gause-Lobatto points where dTy/dz = Ty =0. Since Tjt(z) is 
a polynomial in z of degree fc, ujy is an ^/-polynomial approximation to u(z). 

In order to represent the derivative of the function at the collocation points, Lagrange 


interpolating polynomials of degree N are used, i.e., 

u (*) = J2 **(*)«(**) 

Jk=0 

where the interpolant A*(z) can be expressed in an explicit form: 

\ (~\— (1 - * 2 ) 


(4.14) 


(4.15) 


where cq = cy = 2 and c* = 1 for 0 < k < M. Of course, this Lagrangian polynomial, 
representation provides perfect accuracy for all functions with polynomials of degree < M . 
The p th derivative of ti(z) at the j ih collocation point, in terms of the values of u at all 


the collocation points, is given as 


UT H , t ' _ p 

3-1 i = E a?*"* 


(4.16) 


where u* = u(z*), the Dj k are the elements of the derivative matrix given in Canuto, et 
al. (1987) as 


„ 2>/ 2 + l 

Doo = g — -Vmm 

D jlt = for i + k 

1 c k ( Zj - z k ) 


(4.17) 


= 2(T^J) {oTj = k 
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and, 

O', = (D t „y 

With approximate representations (4.14) and (4.16) for a function u(z) and its derivatives, 
it is possible to convert any ODE into a system of algebraic equations for the [N + 1) 
unknowns ti(zo), ti(zi),.. u(zjy), where zo, zi„. zj/ are the collocation points. These 
(linear) algebraic equations may be solved by several powerful methods [see Wilkinson 
1965], although a principal difficulty is that the matrix associated with these equations is 
not very sparse. 

Clearly, in order to apply a spectral method with the Chebyshev polynomials, the phys- 
ical semi-infinite domain 0 < r < oo must first be transformed onto the finite Chebyshev 
domain, |z| < 1. For boundary layer flows, the two mappings which have been success- 
fully employed in transforming the infinite domain (0, oo) onto the half interval of the 
Chebyshev domain (0, 1) are the exponential mapping (Laurien & Kleiser 1985) and the 
algebraic mapping (Herbert 1984). These mappings, which cluster the collocation points 
near the “wall” at r = 0 and become more widely spaced away from the wall, are inefficient 
in subsonic round jets where steep gradients occur in a very thin critical layer of thickness 
0(i?e”S) located in the neighborhood of the jet radius, not on the jet centerline. The thin 
critical layer in jets requires sufficiently dense collocation points in this layer. In fact, a 
difficult numerical task that arises from the coupled system of ODE is that any numerical 
scheme must be capable of capturing accurately the eigenfunctions and their derivatives 
up to second order, simultaneously. The viscous and conduction terms, which contribute 
to the second derivatives, are included to eliminate any singular behavior of inviscid neu- 
tral eigenmodes. We utilize two different mapping methods for the stability calculation 
of compressible subsonic jets. The first method maps the infinite physical domain onto a 
finite Chebyshev domain. The second approximates the unbounded physical domain by 
a finite domain which is mapped onto the Chebyshev domain. We now outline these two 
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mappings. 


Mapping A 


The first mapping maps from r = 0, oo into z — 0, 1 via 


z — 


h , 1 , r”, 

2FFM + 2 expI “<T l 


(4.18) 


The parameters /i I 2 and n are adjustable constants which control the amount of stretching 
in the r direction. Note that when n = 1 and = 0, the mapping reduces to the classical 
exponential mapping. For a given /1 and I 2 , n controls the number of grid points in the 
neighborhood of the critical point (near r = 1) When the value of n increases, more nodal 
points are clustered wound the critical layer. The parameters I 2 and n redistribute some 
of nodal points near r = 0 to the region where strong gradients exist (near r = 1). For a 
given number of collocation points, the concentration of nodes in the physical domain is 
governed by the metric coefficient: 


z' = 


dz 

dr 


h. _ n r n-ir_r\ 

2 (r H-/x) 2 h P[ / 2 1 


(4.20) 


The higher the value \z*\ in the neighborhood of the critical layer, the greater the concen- 
tration of nodal points. As a caveat, in order not to decrease the exponential convergence 
rates associated with a spectral method, a mapping, r = r(z), which maps the physical 
domain onto the Chebyshev domain must not be too complicated [Canuto et al.]. In other 
words, if we expand 

»■*(*) = X)a t T*(z), 

jfc=0 

the truncation error, e, which arises from using the M + 1 modes to represent r, is 


e — 


00 


k=M 


This error e must be so exponentially small that it will not affect the accuracy of a 
spectral method. In (4.18), the point r = 0 is mapped onto z = 1, while the point r — 00 
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is mapped onto z — 0. The point z = 0 is not used as a collocation point, and we used odd 
polynomials expansion so that the boundary conditions at this point are automatically 
satisfied. 

In general, a function and its higher order derivatives which are regular at infinity in 
the physical domain r might not necessarily be regular in the computational domain z 
(Spalart 1984). In order not to destroy a spectral accuracy, the function and all its higher 
order derivatives at z = 0 must be sufficiently well behaved so that the formula (4.8) for 
the spectral decay of the spectrum a * is valid. We now examine the behavior of u(z) as 
z — ► 0. As r — ¥ oo, z vanishes algebraically according to the first term on the right hand 
side of (4.18). Since a disturbance decays exponentially (2.38) - (2.42), we can conclude 
that all the derivatives of the disturbance — ► 0 as z — * 0, i.e., 

lim — ► 0 for all j (4.19) 

*-* o dz* 

Although no studies have been done here, we suspect that when a function is properly 
resolved, mapping A will preserve the rapid convergence rates associated with a spectral 
method. Further, the algebraic term of this mapping provides an adequate resolution of 
nodal points far away from the origin. 

The cylindrical geometry of a round jet inherits a coordinate singular point at r = 0 
which requires a little care. To remove this singularity, we imposed explicitly the boundary 
conditions at r = 0. These boundary conditions, which are independent of the viscosity, 
can be obtained from (2.38) - (2.42). For f) = 0 or = an even integer, they are 


du 

= dp = 

dp 

dr 

dr 

" dr 


(4.21) 


and, for /? = odd integer 


dv 



(4.22) 


Mapping B 
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For a disturbance which is confined to the critical layer at r — 1 and which vanishes 
rapidly away from the origin, an alternate approach is to transform a truncated finite 

physical domain 0 < r < i?max onto the computational interval 0 < z < 1 according to 

z — cos 8 


y = 1 — 28 n 


r 


f sinh[r(y - B)] ) 
c \ sinh(rB) J 


(4.23) 


where 


B — —In 


[1 + (^ r "" l)( r c/^max) 


, 0 < T < OO, 


2 T Ll + (e T - 1) (r c / i?max) . 
t is the concentration parameter which varies from zero for no stretching to high values 

which concentrate more grid points near r c , and r c is the approximate location of the 
critical layer. Figure (4.3) shows the distribution of collocation points for various values 
of r, [See Anderson et al. 1984 for a part of this mapping]. Since the slope, dz/dr f is 
large in the region of the critical layer, the mappings (4.23) provide high resolution for a 
suitable value of r. Unfortunately, because dz/dr vanishes at Rmax> any disturbance under 
this mapping has zero slope at i?max> unless the asymptotic behavior of this disturbance 
is implemented appropriately (see remarks in section 4.1) into the discretization. As a 
consequence, this mapping will work only for the case where a disturbance is confined 
mainly in a narrrow region and decays rapidly. Extensive sensitivity studies were per- 
formed to determine the effect on the eigenmodes due to the zero slope approximation of 
a disturbance at We found that when aRm** is sufficiently large, typically 4 to 7, 


the stability characteristics are rather insensitive to the value of R max as long as there is 
an adequate clustering of points in the critical layer. This insensitivity of the numerical 
results can be attributed to the sufficiently rapid decay of a disturbance. According to 
the asymptotic solutions (2,38) - (2.42), a disturbance is of order RmJz exp(“ai2max) at 
r == Rmax , which is small enough to be approximated by zero derivative whenever aR max 
is sufficiently large. In this work, r is between 5 to 10, while r c ranges from 0.95 to 1.2. 
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In either mapping, as a result of using either odd or even polynomials over the half 
interval, the spectral derivative operators (4.16) need to be modified accordingly. The 
matrix elements of odd and even derivative operators become 


Of.(odd) = Z>J,(full) - 

(4.24) 

Dj t (even) = D’,(ftlll) + 0£j,_j(full) 

(4.25) 


where 


N = 


M + l 

2 


J = 0,1, 


N 


*: = 0, 1, - - • , iST 


The superscript p denotes the p ih derivative, M is an odd integer, and N is the actual 
number of collocation points in the half Chebyshev domain. The derivatives in the phys- 
ical domain can be expressed in the computational domain by using the chain rule of 


differentiation. For example, 


u(r) = u(z(r)) 


<Pu 

dr 2 


du 

dr 


du dz 
dz dr 


_ <Pu.dz. 2 

= d^ [ d^ } 


+ 


du d 2 z 
dz dr 2 


(4.26) ' 

(4.27) 


where 


dz _ dz dO <fJ7 
dr dO d$ dr 

d 2 z _ dz dy 2 <P9 <Pz dy 2 .d<^. 2 dz d9 cPy 

dr 2 dO dr dy 2 d$ 2 dr dy d& dy dr 2 


(4.28) 

(4.29) 


In an implementation of a spectral collocation method, the spectral derivative operators 
are replaced by full matrices, non-constant coefficient terms are represented by diagonal 
matrices, and the boundary conditions (when exists) are imposed by replacing the ap- 
propriate rows of matrix. The spectral matrix representation of the stability equations 
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converts them into the generalized algebraic problem 

AX = <rBX (4.30) 

where A and B are complex general square matrices of size N*; N* is 5 xN for the primary 
stability problem, 10 xN for the secondary subharmonic analysis, and N is the number of 
collocation points in the computational domain. 

There are generally two different classes of numerical approaches for computing the 
eigenvalues <j of (4.30), global and local methods. Global methods, which require no 
initial guess on the eigenvalues and compute all the eigenvalues of the discretized system, 
are based on the QZ algorithm (Wilkinson 1965). For a matrix of dimension N * , it 
requires 0 (N* 2 ) storage and 0(A"* 3 ) operations. Local methods, however, require a starting 
guess for the eigenvalue and an iteration scheme based on an accelerated inverse Rayleigh 
iteration method (Wilkinson). The global algorithm requires much more computer time 
and storage than the local method since they compute the whole spectrum of eigenvalues 
of the discretized system. Therefore, in order to maximize the numerical efficiency, it 
is desirable to minimize the number of Chebyshev polynomials used in a global method. 
On the other hand, the main shortcoming of a local method is that it converges to the 
true eigenvalue only if an initial guess is sufficiently close to it. As a rule of thumb, in 
the absence of any prior guidance on the eigenvalue, we select the physically desirable 
eigenvalue from a global method and use it as input to a more refined local method 
(using more collocation points) to obtain a more accurate eigenvalue and its corresponding 
eigenfunction. 

For the tanh mean profiles of (2.28) - (2.30), we found that the eigenvalues obtained 
from the global method and the local method agreed to at least six decimal places when 
the same number of grid points were used. The typical rapid convergence of complex phase 
velocity C corresponding to the most unstable eigenvalue of a near maximally amplifying 
disturbance is reflected in the values below: 


84 


N 

c R 

Ci 

15 

0.57644 

0.18860 

21 

0.57656 

0.18866 

31 

0.57656 

0.18866 

41 

0.57656 

0.18866 


The above results, obtained using Mapping B> are for a — 3.5, /? — 1, Mj — 0.1, Re — 
10 6 ,T* = 1 and 9 = 1/15. The corresponding value of C obtained by an independent 
inviscid finite difference code (see Section 4.2) using 200 grid points is C — 0.57656 + 
0.18867t. 

For a near neutral wave the convergence of the unstable eigenvalue with N is observed 
to be much slower than that of a maximally amplified wave. For example, the least stable 
eigenvalue for a near-neutral wave of a heated and compressible jet [a = 5.8, /3 = 0, My = 
0.8, Re = 5000, T* = 2 and 9 — 1/15] is found, using Mapping to be 


N 

Cr 

Ci 

35 

0.4821 

0.0096 

45 

0.4840 

0.0109 

65 

0.4840 

0.0100 


This slow convergence is attributed to the near singular behavior of the eigenmodes in 
the vicinity of a critical point at high Reynolds number flows (see section 2.7). In fact, 
as heating and compressibility increase for high Reynolds numbers flows, we have found 
that the structures of the near-neutral eigenmodes become increasingly more difficult to 
resolve. 

Figures (4.4) - (4.9) exhibit the typical transverse distributions of the shape and phase 
of the most amplified and near neutral eigenfunctions. These eigenfunctions are normalized 
such that the maximum magnitude of the eigenfunction of the streamwise velocity is 0.5. 
The pressure eigenfunction is seen to exhibit the ^smoothest” structure. This behavior is 
sensible since the pressure perturbation has the weakest form of singularity at the critical 
point. Numerically, it is precisely this relatively smooth behavior of the pressure that 
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makes it the most useful variable to work with in connection with a Rayleigh equation 
(when one exists). It is interesting to note that for a near maximally amplified wave, the 
streamwise perturbation u has a typical distinct peak and the radial component t; has a 
plateau- like structure at the critical layer. For a near neutral wave, the peak structure of u 
perturbation changes drastically to a shape with two local maxima and is accompanied by 
a significant change of the phase, while the plateau feature of the v component transforms 
to a smooth shape with a single maximum point. Note also that the amplitude of a highly 
amplified u perturbation is large at the jet lipline and becomes progressively small as the 
wave approaches the neutral point. 


4.2 Some Remarks on Spectral Methods 

1. The imposition of the asymptotic boundary conditions (see section 2.6) forbids the 
use of a global method because the eigenvalue a then enters nonlinearly into the 
eigenvalue problem through the boundary conditions. In addition, since the analytic 
far field behavior of a disturbance in a periodic base flow is still unknown, it is 
suitable to use the set of mappings (4.18) and (4.23) with a sufficient number of grid 
points to resolve the disturbance structure in the critical layer. 

2. The extension to spatial stability analysis can be done without much difficulty. Since 
the basic instability mechanism is essentially inviscid, the viscous terms which con- 
tribute to a quadratic nonlinearity in a, can be neglected in a globed method for 
spatial stability problems. This is because a 2 ~ d 2 jdz 2 and in the critical layer 
(where the effects of viscosity and heat conductivity are small), this derivative is 
negligible compared to d 2 jdy 2 (i.e., the boundary layer assumption). This inviscid 
spatial analysis can be reduced to a generalized eigenvalue problem for a. The phys- 
ically meaningful a, usually the most unstable mode selected from a global scheme, 
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is used as a starting guess for the eigenvalue of a more refined viscous local code. 


4.3 Fourth-Order Finite-Difference Method 


In order to have confidence in the spectral method, the eigensolutions of the viscous 
problem as viscosity vanishes, or equivalently as the Reynolds number becomes infinite, 
were checked by an independent fourth order inviscid finite difference method. Eigensolu- 
tions of the linear stability problem of unidirectional round jets have been obtained using 
finite difference methods by Paragiri (1985) for inviscid analysis and by Jarrah (1989) 
for viscous modes. For completeness, a brief description of the method is discussed here. 
In the absence of viscosity (and conductivity), the governing disturbance equation for an 
inviscid compressible parallel flow can be reduced to a single Rayleigh-type second order 


ODE for the pressure disturbance (2.48) 


d 2 p ri 2 uq 
dr 5 ' :j r uq — C po\ 


Po 






a 2 (u 0 - rypoMj - - a 


P = o. 


(2.48) 


Since tig and p' 0 vanish as r — ► 0 and as r — ► oo, the asymptotic solutions are 


p = C\Ip[ctryJ\ - M?*>o(uo - C) J ) as r — ► 0 (4.31) 

p = CjKp(aryJl - Mjp^uo ~ C) 2 ) as r oo (4.32) 

where Ip and Kp are modified Bessel functions of the first and second kind, respectively, 
and Ci and Ci are constants. Numerically, the asymptotic far field boundary condition is 
satisfied at approximately r = Rm ax = 4.0 and the boundary condition on the jet axis is 
satisfied at the first two grid points where constant mean flow prevails to a good degree 
of approximation. Details and some sample calculations may be found in Paragiri (1985). 
A fourth order central difference scheme is employed which translates the ODE (2.48) 
together with the boundary conditions (4.31) and (4.32) into a pentadiagonal system of 
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complex algebraic equations. The complex eigenvalue C is then found by a local method 
based on the Newton-Raphson scheme where iterations on the eigenvalue are performed 
until a prescribed convergence criterion is met. In this work, the eigenvalues are iterated 
to five decimal places for highly amplified modes, and to three decimal places for near 
neutral modes. Note that a physical eigenvalue can easily be distinguished from a spurious 
(unphysical) one because the latter does not converge as the number of computational grid 
points is increased. 


4.4 Some Remarks on the Linear Inviscid Calculation 


The inviscid pressure equation (2.48) has a regular singular point whenever the 
mean flow velocity equals the phase velocity of a neutral wave. This singularity often 
creates numerical difficulties. For sufficiently amplified modes, the imaginary part of the 
complex phase velocity C removes the singularity from the differential equation along 
the real axis and hence poses no numerical problem. However, for eigenmodes that are 
merely slightly unstable, the almost singular nature of the differential equation requires 
some care. To make matters worse, the secondary instability calculations require an 
accurate determination for all the primary disturbances (not just the pressure) and their 
higher order derivatives. As shown from section (2.7), compressible flows introduce density 
fluctuations whose leading order behavior in the neighborhood of the critical point is 


Pb = 


f (a 2 r» + l 2 r c K + [3*r e 3 - (a 2 r e * + /?*)]«' 1 
1 a 2 r?u' s (r - r c ) J 


H 1- Kp a Iog(r - r c ) (2.56) 


This singularity, which is a pole of order one, may create a change in density fluctuation 
across the critical layer that is too rapid for an accurate determination of an inviscid 
near-neutral eigenmode. Therefore, we conclude that although the instability is predomi- 
nantly inviscid, the inclusion of viscosity is mandatory in secondary instability calculations 
because the primary wave is very nearly neutral. 
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• n — 6 

• n = 3 
° n— 1 



Fig 4.2 (a) and (b) a as a function of the radial variable r for various n. Crosses 
are the collocat, on points. (c) collocation point, on r. The number of collocation 
pomts is 35 using Mapping- A with /, = 7 and I, = 3. 






Fig. 4.3 (a) z as a function of the radial variable r for various stretching parameter t 
(b) collocation points on r. The number of collocation points is 21 using Mapping-B. 
Crosses are the collocation points. Bullets ( t = 5), circles (r = 10), and squares fr 
= !)• 


91 




Fig. 4.4 Distributions of the |eigenfunction| of a near maximally amplified distur- 
bance for Re = 1.0 E + 6,T, = 1.0, Mj = 0.1,0 = 1/15, a = 3.5,/? = 1, and the 
complex phase velocity C = (0.577,0.189). 







phase of w phase of p 

Fig. 4.5 Distributions of the phase of the eigenfunction of a near maximally ampli- 
fied disturbance for Re — 1.0 E + 6 ,T, = 1.0, M, = 0.1,0 = 1/15, a = 3.5,/? = 1, 
and the complex phase velocity C = (0.577,0.189). 












Fig. 4.6 Distributions of the |eigenfunction| of a near maximally amplified distur- 
bance for Re = 1.0 E + 3,T t = 1.0, M, = 0.1,0 = 1/15, a ~ 3.5,/? = 2, and the 
complex phase velocity C = (0.557,0.142). 






Fig. 4.7 Distributions of the phase of the eigenfunction of a near maximally ampli- 
fied disturbance for Re = 1.0 E + 3 = 1.0, M, = 0.1,0 = 1/15, a = 3.5,/? = 2, 
and the complex phase velocity C = (0.557,0.142). 
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Fig. 4.8 Distributions of the |eigenfunction| ot a near neutral disturbance for 

Re - S.OE + 3 \,T. = 1.0, Af, = 0.8 , 6 = 1/15, a = 6.55,/? = 0, and the complex 
phase velocity C = (0.569,0.001). 
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Fig. .1.9 Distributions of the phase of the eigenfunction of a near neutral disturbance 

for Re - 5.00 + 3 ,T. = 1.0, M, = 0.8,8 - 1/15,0 = 6.55,0 = 0, and the complex 
phase velocity C = (0,569,0.001). 
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Chapter 5 

LINEAR STABILITY OF PARALLEL JETS 

The linear stability of a parallel jet is important in its own right. In order to be 
able to manipulate the downstream evolution of jet shear layers, the early stages of jet 
instability must be understood. We reiterate that linear analyses can reveal accurately 
not only the stability characteristics of the initial stages of development of the large-scale 
instability wave, but also, rather surprisingly, it can predict the cross-space structure of the 
perturbation, even when the local disturbance level is as high as 24% of the jet exit velocity 
(Petersen k Samet 1988). This linear analysis, being the first cascade of instability, is 
used as input to fully nonlinear numerical simulations, weakly nonlinear analysis, and 
secondary instability theory. This section considers some features of the linear stability of 
subsonic heated round jets. The different parameters affecting jet instability mechanism 
are discussed separately in order to elucidate their individual influences. The important 
parameters are the jet Mach number Mj y the temperature ratio T*, the jet momentum 
thickness 0 , wavenumbers (a, /?), and the Reynolds number Re . 


5.1 Viscosity Effect 

There are now a number of investigations on the inviscid stability characteristics of 
round jets, although studies on the viscous instability analysis remain scarce. Batchelor 
k Gill (1962) pioneered a theoretical framework on the inviscid stability of a top hat 
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velocity profile which characterizes the mean flow close to the jet exit. This top hat 
velocity, which ignores the presence of the jet nozzle, as is customary in all parallel flow 
analysis, provides a closed form solution of the stability characteristics of a cylindrical 
vortex sheet. By using more realistic profiles, which depend continuously on r, Michalke 
(1971) and Mattingly ic Chang (1974) obtained a good agreement with the experimental 
results of Crow & Champagne (1971). Morris (1976) went further to investigate the effect 
of viscosity on the instability of various jet mean velocity profiles for different values of the 
jet parameters. He found that for jet velocity profiles with thin shear layers, i.e., close to 
the jet nozzle, viscosity has a slight stabilizing effect on the growth rate for flows above a 
Reynolds number of 1000, based on jet centerline velocity and jet nozzle radius. In other 
words, viscosity for instability modes in jets is predominantly passive in character, and 
can be ignored for flows at high Reynolds number (except for neutral waves). 

Typical results exhibiting the influence of Reynolds number on the linear (temporal) 
instability of subsonic heated round jets with the tanh base profiles of (2.28) - (2.30) are 
summarized in Figures (5.1) - (5.6). In particular, the effect of Reynolds number on the 
full spectrum of eigenvalues using 31 Chebyshev polynomials for a = 3.66 and 6 = 1/15 are 
plotted for various jet parameters in Figures (5.1) - (5.3). We caution the reader that the 
discretized eigenvalues in these figures, except for the unstable ones, are only “roughly” 
determined since their associated eigenfunctions did not adequately satisfied the far field 
boundary conditions. Nevertheless, these eigenvalues may provide a rough guide on the 
the full spectrum of the primary instability problem. The salient features that can be 
extracted from these figures are 

1. As the viscosity tends to zero (Re — ► oo), the spectrum consists of a string of neutral 
modes and a complex conjugate pair of eigen values. 

2. The string of eigenvalues, which probably constitutes part of the continuous spec- 
trum, generally becomes increasingly stable with decreasing Reynolds number. The 


c-2. 
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inclusion of viscosity (at high Reynolds number) only slightly stabilizes the unstable 
(discrete) eigenvalue, and breaks down the symmetry of an in viscid complex conju- 
gate pair. In other words, viscosity has a considerable effect on damped waves but 
exerts only a slightly stabilizing effect on unstable (inflectional) waves, 

3. The “continuous modes”, whether viscous or inviscid, have phase velocities lying 
between the minimum (free stream value) and the maximum (jet center line value) 
of the jet mean flow. 

Observation (1) is perfectly consistent with the well-known properties of the reduced in- 
viscid equation (2.48). Since this equation admits a complex conjugate pair of eigenvalues, 
it has either neutral modes or growing modes which have corresponding “damped” modes 
moving with the same phase velocity. Furthermore, since the mean velocity profile used in 
this work has only one generalized inflection point, the spectrum is thus restricted to an 
assembly of a single unstable mode, its conjugate damped mode, and the remaining neu- 
tral modes. The fact that an inviscid equation describes the same instability mechanism 
as the viscous problem in the limit of vanishing viscosity is reinforced by (2). Although 
the continuous spectrum of free shear flows is still not completely understood, observation 
(3) suggests that the complex eigenvalues of the highly damped continuous modes move 
toward the real axis. This observation is again consistent with the theoretical results of 
the continuous spectrum of boundary layers (Grosch ic Salwen, 1978); they found that 
the phase velocity of the highly damped continuous mode matches the mean flow velocity. 
In view of the physical insignificance of the damped waves compared to the rapidly ampli- 
fied unstable waves of round jets, the continuous spectrum will not be pursued further in 
this work. Nevertheless, the continuous spectrum obtained from a numerical discretiza- 
tion scheme can influence the results of an unstable mode. The inclusion of viscosity will 
generally increase the rate of numerical convergence of the unstable mode by moving the 
continuous spectrum away from it. This implies that viscosity will play an important role 
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in the convergence of a near neutral unstable eigenvalue (see Santos, 1987). 

In Figures (5.4) - (5.6), the growth rate and phase velocity as functions of wavenumber 
are plotted for 9 = 1/15. These figures show the typical influence of Reynolds number on 
the unstable mode of subsonic heated round jets. Note that the eigenvalues asymptotically 
approach the inviscid limit at a Reynolds number of about 5000. In contrast, the inviscid 
limit of a planar mixing layer occurs at a Reynolds number of about 200 (Jarrah 1989). 
This means that viscosity plays a more important role in the instability of round jets 
than in mixing layers. This is because in a mixing layer, where there are two external 
streams, the viscous solution (see section 2.6) is completely negligible; but this solution is 
not negligible on the jet axis unless the jet momentum thickness is extremely small. Thus 
the viscous solution will affect the eigenvalue through the boundary condition near the jet 
axis, even when the Reynolds numbers is on the order of several thousand. For jet flows 
at Re < 5000, viscosity not only slightly reduces the range of unstable wavenumbers but 
also the growth rate. Viscosity is also seen to have virtually no effect on the phase velocity 
of unstable disturbances, even though the Reynolds number can be as low as 500. These 
results agreed perfectly with the viscous calculations of Jarrah (1989), and are consistent 
with the numerical results for the spatial and viscous instability calculations of Morris 
(1976). 


5.2 Mach Number Effect 

The new quantity in compressible flows is the density fluctuation, which occurs 
simultaneously with the velocity and pressure fluctuations. The new parameters which 
will generally influence compressible flows are the Mach number My, the temperature ratio 
T*, the Prandtl number Pr, and the ratio of the heat capacity coefficients, 7 = C p /C v . In 
this section, the influence of Match number on the stability of jet shear layers is given, the 
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effect of temperature ratio will be discussed in next session, while the last two parameters 
are assumed to be constants [7 = 1.4, Pr == 0.72], and hence, will not be addressed in 
this work. The Mach number can influence the instability via the base density profiles 
as well as through the disturbance equations. For example, in the unheated jet, the jet 
Mach number appears explicitly in the disturbance equations as well as in the base density 
distribution. 

Although the stability of incompressible shear layers has been studied extensively, a 
systematic determination of the role of Mach number on the instability of jet shear layers 
is still unavailable. In the case of mixing layers, on the other hand, there are now several 
fairly complete studies [see Gropengiesser 1967, Jackson k Grosch 1988]. In his spatial 
stability calculations, Michalke (1971) found that increasing Mach number reduces the 
range of unstable waves as well as stabilizes the flow. Furthermore, he found that the 
phase velocity is sensitive to Mach number only for lower frequencies (or, wavenumbers) 
while in contrast, the growth rates become almost independent of Mach numbers. His 
results also show that the stabilizing effect due to Mach number is slightly more significant 
for an axisymmetric mode than for the first helical mode. 

Eckart (1963) extended the Rayleigh stability criterion and Howard’s (Howard 1961) 
semi-circle theorem to a compressible, adiabatic plane jet. By using a different approach, 
Blumen (1970) found that the influence of aMj on the complex wave velocity of an unstable 
mode is the same as that of the Richardson number in Howard’s results. He concluded 
that compressibility reduces the range of unstable wavenumbers and stabilizes the mixing 
layer. We will show that Howard’s semicircle analysis can similarly bound the complex 
wave speed and illuminate the role of compressibility in jet shear layers. We will do this 
by generalizing Blumen’s approach to a round jet. 

The linear disturbance equations for an inviscid compressible jet can be reduced to a 


single equation for the pressure perturbation p(r): 




2uq Pol d P 


dr + 


\J - »>f2 ? 


a 2 (uo - CypoMj - — a 


p = 0. (2.48) 
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We multiply the above by r/{po( u o — C) 2 } to make it self-adjoint. The equation in the 
self-adjoint form becomes 
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where qi - r [a 2 (u 0 - C) 2 p 0 M 2 - fi 2 /r 2 - a 2 ] / {p 0 (uo - C) 2 } . After multiplying (5.1) by 


the complex conjugate of p and integrating over the range of r, we find that 
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Separating (5.2) into the real and imaginary parts, we obtain 


(5.2) 
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where 


q 2 = 
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C = C R + iCr 


For unstable waves (Cj > 0) which vanish at infinity, (5.3) and (5.4) can be simplified to 

POO POO 

/ uoq 2 dr = Cr q 2 dr (5.5) 

Jo Jo 

POO POO POO 

I u 2 q 2 dr = (Cr + C 2 ) I q 2 dr + a 2 Mf\p\ 2 dr (5.6) 

Jo Jo Jo 

If we let a and b be, respectively, the maximum and minimum of the base velocity [i.e., 
a < u 0 (y) < 6], then (5.5) and (5.6) can be reduced to 
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0 > |[Cn - \{a + 6) 2 ] + Cj - [|(a - 6)] 2 } jH q 2 dr + J~(aM j ) 2 \p\ 2 dr (5.8) 

The above expression can be simplified to the equation governing the Howard’s semicircle 
theorem: 

{[C* - i(« + J)] ! + C?} < [l(a - ft)] 1 (5-9) 

Hence, the stability of inviscid compressible round jets, here a = 0 and 6 = 1, can be 
reduced to exactly the same form as in Howard’s investigation of the stability of variable 
density parallel shear flows. 

Therefore, the complex wave velocity for an unstable mode of jet shear layers is also 
restricted to the semicircle which has the range of uq for diameter, just as in unbounded 
planar mixing layers. Furthermore, since equations (5.5) and (5.6), which govern the 
instability of round jets, are exactly the same as equations (3.4) and (3.5) in Howard’s 
results, we can conclude that compressibility stabilizes jet instability modes. 

For the tanh profile, with 9 = 1/15, the typical role of subsonic Mach number in jet 
instability for the axisymmetric and the first helical modes is illustrated in Figures (5.7) 
and (5.8), respectively. Compressibility is found to be stabilizing and reduces the range of 
unstable wavenumbers. For example, when the jet Mach number, My, is increased from 0 
to 0.8, the maximum growth rate of the axisymmetric and the first helical modes for an 
unheated jet is reduced from 0.68 to 0.55, and from 0.66 to 0.55, respectively. Compress- 
ibility also reduces, but only slightly, the phase velocity of waves whose wavenumbers are 
smaller than the maximally amplified wavenumber, and has no apparent effect on higher 
wavenumbers. In other words, the phase velocities of disturbances, whether axisymmetric 
or helical, are only weakly dependent on the subsonic Mach number. 

In the limiting case as a — ► 0, these figures also show that long waves are unaffected 
by compressibility, i.e., waves with axial symmetry travel with the base flow velocity on 
the center line, whereas helical modes travel with half that velocity. Our calculations 
are consistent with the results of Jarrah (1989) and are in reasonable agreement with 
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the spatial stability calculations of Michalke k Hermann (1982). Physically, according to 
Blumen (1970), the mean flow energy, which provfdes the energy for an unstable mode, is 
diminished by working against the elastic force associated with a compressible medium, 
and consequently the growth rate of an unstable disturbance is reduced. 
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5.3 Heating Effect 


In the previous section, we have discussed the effects of jet Mach numbers on 


the instability. An extension of the Howard semicircle theorem provided some analytical 
understanding for the reduction of the growth rates with increasing compressibility. We 
would now like to gain some understanding of the effects of heating. To do this, we shall 
begin by generalizing the vortex sheet problem of Batchelor <k Gill (1962) to a hot jet. In 
order to obtain simple results, we set Mj = 0 in the final expressions. 

As mentioned in section (5.2), the presence of density fluctuations in a zero Mach 
number flow is caused by a variable base density profile which is a function of 0 and 
T*. Although stability calculations have now been successfully carried out for unheated 
jets, studies on the influence of heating on jet instability characteristics remain scarce. 
Michalke (1971) and Michalke & Hermann (1982) are among the few who investigated the 
stability characteristics of subsonic heated round jets by using linear spatial (parallel flow) 
stability theory. They found that the local growth rate of the spatially growing instability 
wave in a zero Mach number flow (i.e., Mj = 0) increases as the total temperature is 
increased. 

In order to shed light on the effect of heating and to increase the understanding of the 
stability characteristics of a viscous heated round jets with continuous velocity and density 
profiles, the solution for a cylindrical inviscid and heated vortex sheet is investigated here. 
The stability analysis of a vortex sheet has a closed form analytical solution, and has 
been given by Batchelor Gill (1962) for the incompressible case. We now consider the 
corresponding compressible case. 

The base profiles of a compressible vortex sheet are 

uq — U a — constant, po = Pa — constant for r <1 

(5.10) 


uq = 0, Po — Po — constant for r > 1 

On both sides of the vortex sheet, given by r = 1, where irrotational flows persist, the 
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pressure perturbations are 


Pa = Ehl0[aT a r) for r < 1 (5.11) 

p b = D h K 0 (ar) for r > 1 (5.12) 

where 

T a = y/l- {U a - C)*Mjp a 

where Ip and Kp are Bessel functions, and Eh and Dh are arbitrary constants. The corre- 
sponding radial velocity can be obtained from the linearized radial momentum disturbance 
equation. They are 


E k Pff(aT a r) 


V a = - 


(5.13) 

(5.14) 


iap a (U a - C ) 

D H K' fi {ar) 

iJic 

where primes denote differentiation with respect to the argument in parenthesis. Let rj be 
the radial displacement of the vortex sheet which assumes a normal mode representation 
of the form 

r] — A/iexpi[a(a; - Ct) + + c.c. (5.15) 

where Ah is a free constant which denotes the arbitrary amplitude. 

The kinematic boundary condition that the interface I is a material surface (i.e., the 
interface is formed by the same fluid particles at all times) can be expressed as 


^ = 0 for r = 1+ 

Dt 


DI 

Dt 


= 0 for r = 1" 


(5.16) 

(5.17) 


where 
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_ = _ + (v + vo) . v 


and, 


7 = I(x,r,<f>,t) 
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Substituting equations (5.10) - (5.15) into (5.16) - (5.17) and after linearizing the 
disturbance equations, we obtain 


• frr xi \ , E k r fi {a T 0 ) 

ta(U a - C)A k = — jjf ttf - 

Pa(U a - C)ta 


iaCA h = 


Dh£p(a) 

PbCia 


(5.18) 

(5.19) 


Furthermore, the dynamic condition that the pressure is continuous at the interface yields 


Ekfy{aT a ) = D h K p (a) 


(5.20) 


(5.18) - (5.20) provide the the equation for the determination of A*, E * and D*. In order 
to have a nontrivial solution, these equations yield 

- •)’ = T.a)2. (5.21) 

C Pa 

where 

. , T , K^a)P,( T„a) 

The complex wave velocity in equations (5.21) cannot be solved analytically for a nonzero 
Mach number. In order to investigate the individual influence of heating, the solution for 
an incompressible heated jet is examined. In this case, an explicit analytical expression 
for the real and imaginary parts of the complex wave velocity can be obtained as below: 

U a 


C R = 


C7 = 


1 + T*Z/0(a) 
U aS J r.Ip(a) 


(5.22) 


(5.23) 


1 + T t Lp{a) 

where T, = Phi Pa- Since the jet centerline velocity is used in the normalization in this work, 
U a = 1. Thus, the incompressible vortex layer is unstable to all disturbance wavenumbers, 
as in the classical Kelvin-Helmholtz problem. Since very short waves for real flows with a 
finite shear layer thickness are stable, this vortex sheet model is applicable to disturbances 
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whose wavelength is much longer than the shear layer thickness. This will be further 
elucidated in subsequent discussion. Since the wavelength of a disturbance is inversely 
proportional to its wavenumber, long wave approximations are, therefore, characterized 
by small wavenumbers and thus low frequencies. 

Using (5.22) and (5.23), the typical influence of heating on the value of Cj and Cr as 

functions of the axial wavenumbers for /? = 0 and /9 = 1 are plotted in Figures (5.9a) and 

(5.9b), respectively. It is found that, except for long axisymmetric waves, heating reduces 
the phase velocity substantially; for example, the phase velocity of an axisymmetric wave 
with axial wavenumber 7.5 is reduced from 0.53 to 0.36. For waves with a approximately 
less than 1.5, this reduction is more pronounced for the first helical mode than for the 
axisymmetric mode. 

In the limiting case, the values of C for long wave approximation (a — > 0) can be 
obtained from (5.22) - (5.23). They are 

C/e = l, C/ = 0, for 1 = 0 (5.24) 

= ITr? Cl = Trk for/? *° ( 5 - 25 ) 

Therefore, a heated vortex sheet is always stable to long axisymmetric modes, but is 
unstable to all non-axisymmetric modes whose growth rates are reduced by heating. For 
a cold jet (T* = 1), long waves with axial symmetry travel with the speed of the center 
of the jet, whereas all helical modes travel with half that speed (Batchelor & Gill 1962). 
The phase velocity of long axisymmetric waves is independent of heating, whereas that of 
non axisymmetric waves is substantially reduced by heating. 

The stability characteristics of a continuous base velocity profile are in general not 
expressible in closed form, and, therefore, must be determined numerically. The phase 
velocity and growth rate as functions of wavenumbers for the tanh velocity profile with 8 
= 1/15 are displayed in Figure (5.7) for (3 = 0 and in Figure (5.8) for /? =1. Note that 
axisymmetric disturbances are more unstable in heated jets than in cold jets only at small 
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axial wavenumbers, and that heating reduces the range of unstable waves. Except for very 
small axial wavenumbers, heating generally decreases significantly the phase velocity; for 
example, the phase velocity of a near neutral axisymmetric disturbance is decreased from 
0.562 to 0.476. While the phase velocity is rather insensitive to Mach number, it is, except 
for long axisymmetric waves, substantially reduced by heating. 

By using a slightly different velocity profile (profile 1 of Michalke 1984), we found that 
these stability characteristics are virtually indistinguishable from those of profile (2.28) - 
(2.30). This insensitivity of the stability characteristics to the details of the (two slightly 
different mean velocity profiles) convinces us that our results are representative of what 
might be obtained for a wide class of mean flows with reasonable shapes. 

By comparing Figures (5.9 a-b) and Figures (5.10 a-b), the phase velocities obtained 
from a tanh velocity profile with 9 = 1/15 and the vortex sheet model are seen to exhibit 
very similar behavior. In contrast, the quantitative values of the amplification factors from 
these two profiles are considerably different (except at small wavenumbers); this shows the 
true limitation of a vortex sheet theory. 

Further, the long wave results for 9 = 1/15 are consistent with Jarrah’s (1989) calcu- 
lations on mean velocity profiles with a different momentum thickness, namely 9 — 1/5. 
For shear flows, as previously mentioned, the stability behavior of a disturbance which 
has a large wavelength in comparison to a characteristic profile length (e.g., shear layer 
thickness) is unaffected by the form of velocity profiles. The validity of the preceding 
statement needs some clarification, though. For round jets, we found that this statement 
is true for axisymmetric modes but not strictly valid for helical modes. To be precise, the 
phase velocities of long waves (i.e., a -► 0) in a tanh velocity profile with a finite shear 
layer thickness, whether axisymmetric or helical, agreed with that obtained from simple 
vortex models. For the corresponding growth rates, this agreement is found in axisym- 
metric modes but not in helical modes. This disagreement is exacerbated as the helical 


mode number and the shear layer thickness is increased. 
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5.4 Mode Number and Shear Layer Thickness Effects 

Since the evolution of instability modes in an axisymmetric jet is affected by 
the shear layer thickness and the mode number, these effects are concurrently discussed 
in this section. The length scale characterizing the effect of shear layer thickness is the 
dimensionless momentum thickness 9 which, in this work, is normalized with the jet radius. 
There have now been a number of theoretical studies on the stability characteristics of 
incompressible round jets. In their theoretical work, Batchelor As Gill (1962) found that 
a top-hat jet velocity profile (i.e, 9 — ► 0), which represents the mean flow very close to jet 
exit, is unstable to a small disturbance for sill axial and azimuthal wavenumbers. Their 
analysis is supported by Plaschko (1979) who found that the thin jet shear layer (i.e., 
small 9) for a slowly diverging jet is unstable to a large number of discrete azimuthal 
modes. The importance of the parameter 9 in the initial evolution of an axisymmetric jet 
is further explored in detail by Cohen (1986). By using a scaling analysis, he found that 
all azimuthal modes will amplify at identical rates whenever 9 ^ 1 and (A x /A^) 2 <$C 1; 
where A* and A^ are the streamwise and the azimuthal wavelengths. This implies that 
the results obtained from the thin jet shear layer are equivalent to the planar mixing layer 
where the stability characteristics depend only on a single parameter a. 

Although the shear layer near the jet nozzle exit is unstable to a large number of az- 
imuthal modes, only the first few are considered since these are the modes that dominate 
the coherent structures experimentally (Cohen 1986). In fact, the geometrical design of 
most laboratory jets, which consists of a small nozzle emanating from a large chamber 
where an acoustic excitation signal is introduced, cuts off higher azimuthal modes. It is 
for this reason that azimuthal excitation of a jet is usually accomplished by a series of cir- 
cumferentially placed acoustic “drivers” . To be more precise, a linearized monochromatic 
pressure wave traveling inside a plane duct, where for simplicity the mean flow in the duct 
is set to zero, before exiting downstream to the jet nozzle is used for illustration. This 
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wave is of the form 

exp (-*o; e t) cos( ^ ) expt x (5.26) 

* L a oo3 h 

where cv e is the excitation frequency, h is the duct height, is the sound speed, y 
is the cross stream direction, x is the downstream direction, and rid is the duct mode 
number. Equation (5.26) shows that the fundamental mode (n<j = 0) will always propagate 
sinusoidally downstream. However, if the excitation frequency w € is small such that oj e < 
oo/hy then the behavior of the acoustic mode along the x axis is purely exponential, 
and the wave is said to be cut-off. Therefore, at a given frequency, only a finite number 
of modes can propagate out of the jet nozzle exit. In order to excite the instability wave 
whose frequency is on the order of the ratio of the dimensional jet exit speed (Uj) and 
the shear layer thickness (0), the excitation frequency of acoustic waves corresponding to 
most laboratory experiments is often quite low. Consequently, only the fundamental plane 
acoustic wave survives through the nozzle exit. This implies that if the nozzle is round, 
then an axisymmetric low frequency acoustic wave, which accentuates an axisymmetric 
instability wave, will be observed at the nozzle exit. Farther downstream where the rolling- 
up of coherent structure occurs, the axisymmetric and the first azimuthal components have 
been experimently observed to predominate the flows. Therefore, the stability calculations 
are focused mainly on these modes. 

The parameter 9 appears explicitly in the velocity profile (2.28) - (2.30). If the flow 
is compressible, 9 also appears as a parameter in the base density via the coupling with 
the velocity profile. In order to single out the effect of the shear layer thickness, the base 
density profile is taken to be constant (i.e., Mj = 0 and T* = 1). In this work, the tanh 
velocity profiles for two different axial locations are investigated, although the analysis 
can be carried out to include profiles at an arbitrary location, provided that the mean 
profiles are known a priori. The values of momentum thickness which characterize the 
two locations are chosen to be 9 = 1/15 and 9 = 1/5. Crighton iz Gaster (1976) reported 
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that these values of 9 correspond to axial positions 0.5 and 2.7 diameters downstream 
of the nozzle exit in the experimental measurements of Crow k Champagne (1971). Of 
course, the functional relation between 9 and the streamwise direction, which depends on 
Mj and T*, will vary from experiment to experiment. 

The growth rate and phase velocity of the first three azimuthal disturbances for a cold 
and incompressible jet for 9 — 1/15 and 9 = 1/5 are summarized in Figure (5.11) and 
Figure (5.12), respectively. In Figure (5.11), the jet shear layer is dispersive to waves 
longer than the most amplified wave and relatively nondispersive for shorter waves; ft = 1 
produces the largest range of (approximately) nondispersive waves. The phase velocity of 
axisymmetric modes decreases with a from unity to 0.56. In contrast, the phase velocities 
of helical modes (ft = 1, 2) increases with a from 0.5 to 0.56. In particular, the phase 
velocity for a neutral wave is approximately 56 % of the jet centerline velocity, and it 
is independent of the mode numbers when the shear layer is very thin. An important 
consequence of the almost nondispersive character of azimuthal modes over a broad range 
of wavenumbers is that it allows waves to interact resonantly. We shall show in Chapter 6 
that this resonance wave interaction phenomenon is, in fact, a mechanism for the strong 
growth of a subharmonic disturbance. Further, the axisymmetric (ft = 0) mode has the 
highest maximum growth rate of 0.680, the ft = 1 mode amplifies at a slightly lower ^ate 
of 0.663, while the most unstable rate of ft = 2 is 0.595. This is reasonable, since as 9 — ► 
0, the results of the tanh profiles tend to that obtained from vortex sheet analysis where 
all azimuthal modes amplify at the same rate at large enough axial wavenumbers. 

It is interesting to recall that the phase velocity for all azimuthal waves for a simple 
vortex sheet (i.e. for 9 = 0) is at approximately 50% of the jet centerline velocity. This 
shows that although the vortex layer under-estimates the phase velocity of the tanh veloc- 
ity profiles at 9 = 1/15 by about 11 %, it correctly predicts that all azimuthal disturbances 
with high axial wavenumbers propagate at the same phase velocity. 
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Farther downstream at 9 = 1/5, Figure (5.12a) shows the ft = 1, /? = 2, and ft = 0, 
are respectively the most, the second, and the third unstable modes. The maximum 
growth rates of the disturbances with ft = 1, ft =0, and ft — 2 are 0.173, 0.134 and 
0.039, respectively. Therefore, if both the axisymmetric and helical modes have the same 
initial amplitudes, the helical mode will eventually dominate downstream of the end of the 
potential core. The emergence of ft = 1 mode as the dominant instability mode towards 
the end of the potential core is also observed by Michalke ic Hermann (1982), Cohen 
(1986) and Miles & Raman (1988). Batchelor <k Gill (1962) theoretically demonstrated 
that the a fully developed” jet, after the end of a potential core, is unstable only to the 
ft = 1 mode. For the tanh profiles, the phase velocities of axisymmetric modes are seen 
to decrease with wavenumbers from unity to 0.71. In contrast, the phase velocities of 
modes ft = 1 and ft = 2 increase with wavenumbers from 0.5 to 0.67 and from 0.5 to 0.56, 
respectively. Owing to a greater value of 6 , the near-neutral azimuthal modes now travel 
with different phase velocities. 

By comparing Figure (5.11) and Figure (5.12), we found that as the shear layer thick- 
ness grows^not only does the wavenumber of the maximum local growth rate and the range 
of unstable waves becomes smaller, but the range of nondispersive waves is significantly 
reduced. It should be noted that, as pointed out by Petersen (1978), each coalescence 
of coherent structures will result in a bigger value of 0, making the layer more dispersive 
to unstable waves . These results are consistent with the numerical results of Petersen 
(1978), Michalke & Hermann (1982) and Cohen (1986) which are based on spatial stability 


calculations. 
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Fig. 5.1 The effect of Re on the spectrum of linear, parallel flow, temporal eigen- 
values using 31 collocation points for 0 = 1/15, M j = 0, T. = 1, a = 3.66, and 
0 = 0. 
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Fig. 5.2 The effect of Re on the spectrum of linear, parallel flow, temporal eigen- 
values using 31 collocation points for 0 = 1/15, M, = 0, T, = 2, a = 3 66 and 
0 = 1 . 
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Fig. 5.3 The effect of Re on the spectrum of linear, parallel flow, temporal eigen- 
values using 31 collocation points for $ = 1/15, A/, = 0.8, T* = 1, a = 3.66 and 
0 = 0 . 
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Fig. 5.4 (a) Growth rate and (b) phase velocity of unheated, axisymmetric modes 
for various Re\ Re = oo (solid lines), Re = 5000 (dot), Re =500 (dash). Calculations 
are for 0 = 1/15, My = 0, T. = 1. 
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Fig. 5.5 (a) Growth rate and (b) phase velocity of heated, axisymmetric modes for 
various Re: Re = oo (solid lines), Re = 5000 (dot), Re =500 (dash). Calculations 
are for 6 = 1/15, My = 0, T, = 2. 
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Wave Number , a 

Fig. 5.7 Inviscid (a) growth rate and (b) phase velocity of the azimuthal modes 
(3 = 0 and 9 = 1/15 for various Mj and Ty. Mj = 0, T, = 1 (solid lines). 
Mj = 0.8, T t = 1 (chaindot). Mj = 0, T, = 2 (dash). 
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First Azimuthal Mode 




Wave Number, « 

Fig. 5.8 Inviscid (a) growth rate and (b) phase velocity of the azimuthal modes 
P = 1 and 6 = 1/15 for various M ; and Tj. M, = 0, T t = 1 (solid lines). 
Mj = 0.8, T, = 1 (dot). Mj = 0, T, - 2 (dash). 
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Fig. 5.9 (a) Amplification rate (C/) and (b) phase velocity ( C R ) of a vortex sheet 
for 0 = 0 and 0 = 1 modes. f3 = 0, T, = 1 (solid lines). 0 — 0, T t = 2 (dot). 
0=1, T. = 1 (dash). 0 = 1, T t - 2 (chaindash). = M a =0. 
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Fig. 5.10 Inviscid (a) growth rate and (b) phase velocity of the azimuthal modes 
(3 = 0 and 0 = 1 for 0 = 1/15 and M ) = 0. 0 = 0, T t = 1 (solid lines). /? = 0, T, - 2 
(dot). 0 = 1, T t = 1 (dash). 0=1, T, = 2 (chaindash). My = M a =0. 
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Variation of Azimuthal Mode Number 




Fig. 5.12 Inviscid (a) growth rate and (b) phase velocity of the azimuthal modes 
0 = 0 (solid lines), 0 = 1 (dot), and 0 = 2 (dash). Calculations are for 
0 = 1/5, Mj = 0, T t = 1. 
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Chapter 6 

SUBHARMONIC INSTABILITY MODES 


As discussed in Chapter 3, a principal effect of excitation is to generate a primary 
wave which grows in amplitude, and then reaches a finite amplitude (i.e., saturation) as it 
propagates downstream. In turbulent flows, this primary wave is usually called a coherent 
structure. An objective of this work is to investigate, via a secondary instability theory, 
the conditions under which this primary wave can destabilize a subharmonic mode. Since 
subharmonic disturbances generate vortex pairing (and vice-versa), in a physical sense, 
the present chapter focuses on the conditions, as a function of Mach number, temperature 
ratio etc, which enhance, or at least influence, vortex pairing. This process of merging 
(or pairing or amalgamation), is generally believed to be ah important mechanism in the 
spreading of jet flows. 

Although secondary instability theory in round jets is presently still new, its predictive 
capabilities are, for the most part, undisputed in incompressible planar mixing layers and 
boundary layer flows. We now digress momentarily to review some of these known results, 
and outline some distinct features of round jets that are absent in mixing layers and 
boundary layers. 

By using a secondary instability analysis, Kelly (1967) argued that the sequential merg- 
ing of vortices in planar mixing layers can be attributed to the subharmonic instability 
modes in a periodic base flow. These modes, which absorb energy from the mean flow via 
the catalytic role of a periodic fundamental mode, predominate in the flow at approxi- 
mately the location where vortex pairing occurs. Ho & Huang (1982) also demonstrated 
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that after each pairing of coherent structures, the peak frequency of the energy spectrum 
was halved, and the merging process was found to be at the location where the sub- 
harmonic modes become the dominant instability modes. Their findings reinforced the 
predictive capability of a secondary instability theory. 

Further, Herbert (1983, 1984, 1985) demonstrated that a secondary instability analysis 
can unravel some of the fundamental routes to transition in incompressible wall-bounded 
flows. In a particular route, he has shown that the essential flow physics associated 
with the formation of staggered lambda vortices in boundary layers, although seemingly 
formidable, can be revealed by a linear secondary instability theory. It may be remarked 
that boundary layer flows are fundamentally different from shear flows. This difference 
is marked by the presence of a solid wall which constitutes an additional viscous no- 
slip condition on the velocity perturbations. Consequently, the growth of an instability 
mode is induced by viscosity, and hence, the growth rate is small. This means that when 
the viscosity is turned off in wall-bounded boundary layer flows, there are no instability 
modes. In contrast, the instability modes in shear flows are predominantly inviscid, and 
the growth rates are generally large; at least in the sense that the amplitude of the 
fundamental changes by a significant amount in one period or wavelength. 

Jets and planar mixing layers belong to the same class of flows, which are characterized 
by the existence of a generalized inflection point in the parallel flow profiles. However, 
the former have a few distinct features that are absent in the latter. The qualities that 
discriminate jets from mixing layers are: 

1, The radius of a jet nozzle constitutes a length scale, which is used in our normaliza- 
tion, in addition to the shear layer thickness. Thus the shear layer thickness becomes 
a bona-fide parameter for a jet; this quantity can be scaled out for a mixing layer. 

2. The vortex stretching term, which represents the extension or contraction of vorticity 
lines by the strain rate of the flow, is nonzero for azimuthal disturbances. 
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Because of (1) above, the symmetry of the base vorticity about the critical point of a jet is 
broken, and, as a result, the waves are dispersive so that waves with different wavenumbers 
are convected with the flow at different speeds. Furthermore, the criterion for the existence 
of an inflectional mode is now a function of the mean flow and the mode number, /?, of 
the instability mode (Batchelor As Gill 1962). In contrast, the necessary and sufficient 
condition for the existence of an instability mode in a planar mixing layer depends on the 
base flow alone. 

We will employ a generalization of the current incompressible secondary instability 
theories to gain a better understanding of the mechanisms governing the amalgamation 
of coherent structures in compressible round jets. High speed jets, whether for heating or 
for nonzero Mach number, introduce density perturbations. These result in quadratic and 
cubic interactions between the density perturbations and other types of perturbations. 
This cubic interaction results in considerable algebra in the secondary instability analysis, 
and therefore necessitate some simplifications which we will now discuss. 

Recall that a deficiency of linear stability analysis is that it cannot determine a unique 
amplitude A. In this work, the amplitude A is a small arbitrary prescribed quantity 
which corresponds to the maximum amplitude of the streamwise velocity perturbation. A 
consequence of this normalization is that if all other types of perturbations are at most 
the same order as A, then the quadratic interaction between Fi and F 2 will play a more 
important role than the cubic interaction between Fj and itself and F$. This is because, 
for small A, the quadratic interaction is of order As while the cubic interaction is of or- 
der A 2 e . There are, however, possible situations where this definition of A might cause 
inconsistencies with our analysis. These happen in strongly heated jets where the parallel 
stability analysis provides a density eigenfunction that is much greater than the stream- 
wise velocity eigenfunction. As a result, the magnitude of the density perturbation (A/?i) 
will no longer be small compared to the unperturbed density (/>o)- This will invalidate the 
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present stability analysis where the amplitudes of all types of perturbations are assumed 
to be small compared to the mean flow. Nevertheless, these possibilities can be avoided 
by normalizing A so that it corresponds to the maximum amplitude of the density per- 
turbation. This implies that when A is small, the amplitude of all types of perturbations 
will be small too. Therefore, in this work, we shall focus on the cases where A 2 and higher 
order terms can be neglected. The preceding remark is consistent with our representation 
of the flow arising from the finite-amplitude growth of the primary disturbance associated 
with the parallel base flow. Formally, this flow [say, Xi ] may be expanded as 

&(x,r, *) = u 0 (r) + A 2 Uq H HA ju^r) + A 2 u\(r) H J exp(«aix) 

+ |a 2 U 2 + * * * J exp(2iai) H 1 -c.c (6.1) 

where the superscript 1 denotes a correction due to the nonlinear self-interactions of the 
primary wave, and the finite amplitude A can, in principle, be obtained from the so- 
called Landau’s equation (Jarrah 1989). Since the periodic base flow we have modeled 
(i.e., uo + (Aui exp(iaix) + c.c.)) is formally valid to O(A), on the basis of the expansion 
(6.1), any quantity involving A 2 and higher order terms in the secondary stability analysis 
should be truncated. Further, as noted by Kelly, we also assume that the terms of O(A) 
and higher in (6.1) do not eliminate the generalized inflection point associated with the 
parallel mean flow, as otherwise, the perturbation analysis will be invalid. 

Fully developed coherent structures (or roughly, vortices) can be viewed as manifesta- 
tions of a primary instability wave that propagates and grows in the downstream direction. 
Although the quantitative experimental data on the vortex-pairing scenario in round jets 
is presently still very scant, the evidence of the pairing phenomenon can be seen in the 
experimental work of Wille (1963) and Reynolds (1988). In our study , the merging of 
coherent structures is idealized as being composed of two externally excited disturbances: 
an axisymmetric fundamental wave and its subharmonic. The fundamental wave is chosen 
to be axisymmetric for two reasons. First, the axisymmetric geometry of a jet nozzle tends 
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to accentuate an axisymmetric wave near the nozzle. Secondly, the resonance condition 
on mode numbers, which has been established in Chapter 3, can be satisfied exactly for 
a subharmonic mode with an arbitrary mode number. Further, the merging of two ring- 
like coherent structures is idealized as the resonant interaction between two axisymmetric 
waves. And, any appearance of helical motions can be attributed to a helical subharmonic 
resonanting with its axisymmetric fundamental. This idealized model, which assumes 
that other frequencies and fine scale motions play a minor role in the overall flow physics 
describing the merging of coherent structures, provides a blueprint for future studies and 
comparisons. 

To recap, the assumptions used in the secondary instability analysis are: 

1. The magnitude of the subharmonic disturbance is sufficiently small to allow lineariza- 
tion about the primary wave and the parallel base flow. 

2. The mean flow profiles are kept unchanged and the mode shape (i.e., the cross-space 
structure) of the primary wave is given by a linear theory even though the amplitude 
of this wave may be quite large (say, A = 0.10). Some justification of this assumption 
lies in the work of Petersen k Samet (1988). 

3. We shall assume that A is a prescribed constant, although in experiments we find that 
during and after saturation, while the subharmonic is growing, A is a slow function 
of x. In our temporal theory, this translates into A = A(t); we ignore this slow 
dependence on time. 

A consequence of these assumptions is that the stability characteristics of the fundamental 
mode can be approximated by the linear modes of parallel flow analysis. It may be 
reiterated that the disturbances are normalized such that A is the amplitude corresponding 
to the maximum of the streamwise velocity perturbation of the primary wave, and we 
consider a disturbance to be an infinite wave train that is spatially periodic but temporally 
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growing. 


6.1 Resonance Conditions 

In the method of multiple scales, we considered two small amplitude waves of 
arbitrary wavenumbers and growth rates evolving temporally in a parallel shear flow. The 
lowest order nonlinear interaction can occur only through quadratic wave interactions. 
In a perturbation scheme, the flow field at this order may be thought of as the response 
of a linear oscillator to a small forcing function comprised of the quadratic wave (or, 
oscillation) interaction terms. The response of the oscillator to this small forcing function 
is small unless there is a resonant wave interaction. This happens when the growth rate, 
mode number and wavenumber of a forced wave, generated by the forcing term, are the 
same as that of an instability wave of the parallel flow. A consequence of the quadratic 
wave interactions is an enhancement of the transfer of energy from the mean shear flow 
to the forced wave, which after a sufficient period of time, grows to an amplitude that is 
comparable in size to a primary mode associated with the roll-up of a coherent structure. 
Recall that in order to have resonant wave interactions, the two unstable wave modes 
must fulfill the following resonance conditions: 


|ai|Cj(|ai|,|/9i|) + |aj|C/(|a2|,|/9j|) = |c*2 + ai|C7(|<*2 + ai[,|/?i +&|) (3.16) 

a 2 Crt (|ai|,|/9i|) + a^C^da^l, |/%|) = (<*2 + <*i)Crt(|<*2 + ai|, |£i + /3 2 |) (3. 17) 

where the a’s and /?’s can be positive or negative, and are so far arbitrary. These condi- 
tions, upon close examination, can be satisfied by selecting all of the following restrictions: 


(6.1) The wavelength of the subharmonic disturbance is twice that of the primary 
(or, fundamental) wave [i.e, a\ = - 2 a 2 ]. 
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(6.2) The sum of the azimuthal mode numbers of the fundamental and its subhar- 
monic is equal to the azimuthal mode number of the subharmonic alone. 

(6.3) The phase velocity of the fundamental equals that of the subharmonic. 

(6.4) The sum of the growth rates of the fundamental and the subharmonic is the 
same as the growth rate of the subharmonic alone. 


The condition (6.2) can be satisfied exactly by taking the fundamental to be axisym- 
metric, while the subharmonic can be axisymmetric or helical. The condition (6.3) requires 
the waves to be nondispersive and this can only occur when the jet momentum thickness 
is small (i.e., near the jet nozzle). In round jets, where the waves §£. S. dispersive, this con- 
dition is met only approximately. In contrast, the phase velocities in planar mixing layers 
are independent of the wavenumbers and hence guarantee this condition exactly. For the 
tanh velocity profile, the matching of phase velocity can be approximately satisfied in two 
different ways. The first way represents a near neutral fundamental mode interacting with 
a near maximally amplified subharmonic, while the second way involves a near maximally 
amplified fundamental interacting with a subharmonic whose growth rate lies between 
zero for a = 0 and the most amplified value. For clarity, we called these two ways route I 
and route II to subharmonic instabilities, respectively. Route I is physically more mean- 
ingful since a fundamental wave attains a finite- amplitude saturation when it is within the 
vicinity of a neutral point. As revealed in Figure (6.1), an axisymmetric subharmonic can 
interact with the fundamental via route I, while the matching of phase velocities between 
a helical subharmonic and the fundamental can occur through either route I or route 77; 
note that the fundamental is axisymmetric. In Figure (6.1), the subscripts / and s denote 
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the phase velocities of the fundamental and the subharmonic, respectively; symbolically, 
route 7 is represented by an interaction between If and while route II is represented 
by an interaction between II f and II 8 . 

The condition (6.4) can approximately be fulfilled by choosing an almost maximally 
amplified subharmonic and a fundamental wave that is near neutral. Strictly speaking, 
this condition will preclude route II, although sometimes this condition is abandoned 
altogether under the assumption that growth rates are generally “small”. It is worth 
noting that the wave resonant interactions in the bulk of Kelly’s work (1967) on planar 
mixing layers in achieved via route II A consequence of this choice is that the condition 
on the growth rate is not met. The effect on the secondary subharmonic modes for not 
imposing the condition (6.4) will be examined in round jets in subsequent sections. 

In summary, in order for all the resonance conditions to be satisfied simultaneously, 
the fundamental must be a near neutral axisymmetric mode, while the subharmonic can 
have an arbitrary azimuthal mode number and is nearly maximally amplified. 


6.2 Comparison of Approaches 

In view of the absence of any previous studies on the secondary instabilities of 
subsonic heated round jets, two distinct yet complementary approaches are used. One of 
these is a perturbation method and the other is a generalization of normal mode approach 
to a periodic base flow. Typical results from the two approaches for 6 = 1/15 and various 
jet parameters are compared in this section. In Figures (6.2) - (6.4), the solid lines 
represent the method of multiple scales, while the open circles indicate the normal mode 
analysis. In the normal mode analysis, the [relative frequency) = \&u\ — the absolute value 
of the imaginary part of <r 2 defined in (3.45) — refers to the frequency of the subharmonic in 
a frame of reference that moves with the phase velocity of the fundamental wave. When 
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this [relative frequency | is zero, the relative phase velocity in the streamwise direction 
between the fundamental and its subharmonic becomes zero. 

In Figures (6.2 a-b), the growth rate of the secondary wave (i.e., the real part of &2 = 
# 2 r) and its [relative frequency], as a function of the amplitude, A , of a near maximally 
amplified fundamental mode are plotted. These figures represent the case where the phase 
velocity approximately satisfies the resonance condition (6.3), whereas the growth rates 
do not adequately fulfill condition (6.4). A couple of interesting observations can be 
extracted from these figures. First, the normal mode analysis produces a “dominant” 
subharmonic mode whose growth rate increases at a rate proportional to the amplitude 
A of the fundamental and a second mode which suffers a decay as A increases. In the 
perturbation method, we focus on the subharmonic whose growth rate increases with A. 
The growth rate of the dominant mode is increased by about 3% when A changes from 0 
to 0.02. At a larger value of A, the decaying mode is physically insignificant relative to 
the dominant mode and, it will not be pursued here. Secondly, when A ranges from 0 to 
0.02, the growth rate of the dominant subharmonic mode agrees well with that obtained 
from the method of multiple scales. Finally, the phase velocity of the subharmonic wave 
remains synchronized with its fundamental when A becomes finite. Note that in all of 
these calculations, A is still very small; because of the presence of the (periodic) primary 
wave, the subharmonic possesses a growth rate which is larger than what it would have 
in the parallel base flow. 

The previous results are extended to larger amplitudes, A, in Figure (6.3). In comparing 
the two approaches, it can be seen that when A is less than 0.03, they are in agreement. 
However, when A increases from 0.03 to 0.1, the perturbation method overestimates the 
growth rate. We conjecture that this disagreement arises from the inaccuracy of the 
low order perturbation method for large A, which is presumably caused by the limited 
accuracy of the truncated asymptotic expansion (3.21) - (3.24). An obvious suggestion is 



136 


to introduce a second slow time scale Ti — 8*t in order to absorb the parabolic dependence 
of the growth rate on A for large A . Rather than extending the analytical approach in this 
way, we chose to explore the normal mode method whose range of validity is somewhat 
larger than that of the asymptotic method and which also provides an independent check. 
Once the primary wave attains a finite amplitude, the secondary instability of the jet, 
as manifested by the growth rate of the subharmonic, increases. We believe that this 
instability will result in the pairing of the vortex rings associated with the primary wave. 

Figures (6.4 a-b) show the growth rate and the |relative frequency) as a function of the 
amplitude of a near neutral fundamental mode for a heated jet. It can be observed that 
when A increases from 0 to 0.02, the growth rate of the subharmonic mode is augmented 
by about 7%. Further, for this range of A, the two approaches are in good agreement, 
and the phase velocity of the subharmonic mode matches that of the fundamental. 

It is worth noting that the periodic flow destabilizes the subharmonic mode of Figure 
(6.4) more dramatically than that of Figure (6,2). This can be attributed to the fact that 
the results in Figure (6.2) do not adequately fulfill the growth rate restriction (6.4). These 
results are rather typical of our calculations, and therefore, we can conclude that in order 
to achieve optimum enhancement of a subharmonic instability mode, its wavenumber, 
phase velocity and growth rate must satisfy the resonance conditions simultaneously. 

Additionally, in subsequent sections of this chapter, it will be found that in order to 
compare these approaches, it is essential that the phase velocity of the primary subhar- 
monic disturbance matches (or, almost matches) the fundamental. This is because when 
the linear waves are nondispersive, there is no mechanism in the low order perturbation 
analysis to guarantee the destabilization of a subharmonic mode by the periodic flow by 
means of phase-locking with the fundamental. 
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6.3 Compressibility Effect 

In this section, the influence of Mj on the secondary subharmonic instability in an 
unheated compressible jet is presented. It should be remembered that as the jet Mach 
number changes, the shape of the primary mode also changes. Therefore, both this change 
in the primary waves as well the explicit dependence of the secondary wave on Mj will 
alter the instabilities of the subharmonic. [For reference, compressibility stabilizes the 
primary eigenmodes]. The axial wavenumbers of the subharmonic and the fundamental 
are chosen to correspond to a near maximally amplified wavenumber and a near neutral 
wavenumber, respectively. In accordance with the discussion in the previous section, only 
the “dominant” subharmonic is shown. 

A typical influence of Mj on the growth rate and the jrelative frequency | of subharmonic 
instability modes (^2 = 0) and (/?2 = 1) in the presence of a finite amplitude primary mode 
is exhibited in Figure (6.5). With reference to Figure (6.5b), an interesting feature is that 
when the waves are (initially at A =0) slightly dispersive, a finite threshold (fundamental) 
amplitude A^ is required for the matching of phase velocities to be established. After this 
matching is established, the dominant subharmonic mode always travels with the speed 
of the fundamental, and the growth rate of the subharmonic mode increases substantially 
as A increases (see Figure 6.5a). This threshold amplitude generally decreases as the 
difference (at A = 0) between the phase velocities of the fundamental and subharmonic 
decrease. This observation is important because it implies that for small amplitudes of A 
which are less than the threshold value, the low order perturbation series (3.21), which 
assumes that the matching of phase velocities is independent of A, might not be able to 
correctly predict the growth rate of the secondary disturbance when the primary waves 
are dispersive. This is sensible, on physical grounds, since the matching of phase velocity 
means that the two waves travel in unison. Therefore, there is a possibility for the effective 
transfer of energy from the mean flow to the subharmonic mode via the catalytic role of 
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the finite-amplitude primary wave; this is because there is no “phase cancellation” by the 
motion of the subharmonic relative to the fundamental. 

In order to have a representative idea of the effect of compressibility, it is sensible to use 
the total increment of the growth rate of a subharmonic mode from parallel flow analysis 
due to the effect of the primary wave. Let K represent this increment when A changes 
from 0 to 0.1. Conveniently, some of the more important results pertaining to Figure (6.5) 
can be tabulated as below: 


& 

M,- 

A-th 

N 

0 

0.1 

0.011 

0.236 

0 

0.8 

0.004 

0.204 

1 

0.1 

0.006 

0.205 

1 

0.8 

0.011 

0.180 


The above calculations are for Re = 5000, 9 = 1/15, T* = 1, and f}\ = 0. The value of A t h 
is only an estimated one because we did not have enough runs in the neighborhood of A^. 
For an incompressible jet, when the amplitude A changes from 0 to about 0.1, the growth 
rates of the /?2 = 0 and /?2 = 1 are augmented by 0.236 and 0.205, respectively. This is 
about a 40 % increase over their values at A = 0 and suggests that subharmonic instabilities 
play a major role in the development of a jet flow. Furthermore, when A is greater 
than the threshold value, compressibility has a moderate stabilizing effect on secondary 
subharmonic instability modes so that vortex pairing, which results from subharmonic 
instability, is expected to be less important in compressible flows. For example, for the /?2 
= 0, the increment, N, is reduced from 0.236 to 0.204 due to compressibility. 
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6.4 Heating Effect 

In this section, the influence of T t on the secondary subharmonic growth rates in a 
heated (T. > 1) and low Mach number (My = 0) jet is discussed. The comparison between 
heated and unheated flows is again rather subjective since the primary eigenmodes also 
depend on T*. As noted in Chapter 5, heating is found to have a slight stabilizing effect 
on primary modes with axial wavenumbers larger than the most amplified wavenumber. 
As before, the axial wavenumbers of the subharmonic and the fundamental are chosen to 
correspond to a near maximally amplified wavenumber and a near neutral wavenumber, 
respectively (i.e., via Route I). 

Figure (6.6) shows the influence of T. on the growth rate and the |relative frequency! 
of subharmonic instability modes (ft = 0 and ft = 1) as a function of the amplitude 
A of the fundamental mode. In this figure, 9 = 1/15, My = 0.1, and ft = 0 are used. 
Again, for clarity, N refers to the increment in the growth rate of the subharmonic mode 
when A changes from 0 to 0.1, and A t h represents an approximate threshold amplitude for 
the onset of a secondary subharmonic instability mode. The more significant information 
pertaining to Figure (6.6) can be tabulated as below: 

ft T. A th N 

0 1 0.011 0.236 

0 2 0.014 0.206 

1 1 0.006 0.205 

1 2 0.001 0.122 

Severed interesting features can be seen in this figure, using the above table as an 
aid. First, the growth rate of subharmonic instability modes in a heated jet, again, 
increases as A increases. Secondly, while heating slightly suppresses the growth rate of 
an axisymmetric secondary mode (by about 15%), it decreases the growth rate of = 1 
mode considerably (by about 40%). This leads us to believe that heating tends to inhibit 
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vortex-pairing in round jets, and this inhibition is substantially more pronounced for the 
first azimuthal subharmonic mode with ft = 1. Finally, heating increases the threshold 
amplitude required for the onset of a secondary subharmonic with j3 = 0, but it reduces 
the corresponding threshold amplitude for /3 = 1 mode. To conclude, vortex pairing in a 
heated, incompressible jet is likely to be an axisymmetric phenomenon. 


6.5 Shear Layer Thickness Effect 

We recall that jet instability is caused by the induction of the vorticity, which 
is contained within a narrow region in jet shear layers. The physical mechanism for 
an amplified disturbance can be explained in terms of a vortex sheet in a homogeneous 
incompressible fluid (Drazin <k Reid). According to the Biot-Savart law, each vortex in an 
inviscid fluid is carried with the fluid and induces a circulating flow proportional to the 
strength of the vortex line. Specifically, the velocity induced at any point by the vortices 
will amplify the sinusoidal displacement of a vortex sheet. 

In a given parallel velocity profile, the mean vorticity is inversely proportional to the 
thickness of the shear layer. This vorticity of the steady mean flow can be interpreted as a 
source term from which energy is transferred to the unsteady perturbation vorticity. This 
point can be illustrated by considering a two-dimensional incompressible inviscid fluid 
with a steady mean vorticity fio and an unsteady small perturbation vorticity Q f . The 
equations governing the dynamics of the total vorticity and the perturbation vorticity are, 
respectively, 


^(fi» + n') = o 

. d d . , dClo 
( _ + U0 _)n' = - v — 


( 6 . 2 ) 


dr 


(6.3) 
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where 


D 

Dt 



V, 


and uo and t; sire the mean flow velocity and the radial perturbation velocity, respectively. 
These equations reveal that while the total vorticity of a material element carried with 
the fluid is conserved, the perturbation vorticity, which is being convected by the mean 
flow uo, changes along the mean flow because of the “source term” — vdClo/dr. 

Although in secondary instability analysis, the presence of a finite-amplitude funda- 
mental mode will restructure the vorticity distribution, the growth rate of a subharmonic 
mode is, nevertheless, expected to increase as 9 becomes smaller. This growth is revealed 
in Figure (6.7), which is plotted for Re = 5000, My = 0.1, T* = 1, fti = 0 and two different 
values of* 9 = 1/30 and 9 = 1/15. It may be noted that the flow is taken to be unheated 
(T* = 1 ) and almost incompressible (My = 0.1) in order to isolate the effect of the param- 
eter 9 , which appears explicitly in the velocity profile. It is shown that, when A changes 
from 0 to 0.1, the total increment, N, of the growth rates of the axisymmetric subharmonic 
modes for 9 = 1/30 are approximately twice that for 9 = 1/15. These results indicate that 
although the presence of a primary wave redistributes the vorticity of the parallel flow, the 
growth rate of a secondary subharmonic mode is seen to increase proportionally to 1/0, 
just as in the primary stability analysis. This implies that if vortex pairing is a resonant 
interaction between these structures and a small amplitude subharmonic mode, the for- 
mer exerts the greatest influence in reinforcing the latter when the shear layer thickness is 
small (i.e., near the jet nozzle). This means that each successive coalescence of a coherent 
vortex structures not only makes the primary wave more dispersive because of the larger 
shear layer thickness, but also makes it less conducive to pairing. 

In Figure (6.7), it may also be noted that the solid lines represent the case where the 
resonance condition on the growth rate (6.4) is not met, and subsequently, the increase 
in the subharmonic growth rate is less pronounced than in the case where the condition 
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(6.4) is fulfilled. 


6.6 Mode Number Effect 

In this section, we will discuss the influence of the finite amplitude primary ax- 
isymmetric mode on the growth rate and the (relative frequency! of a subharmonic mode 
for various mode numbers. In Chapter 5, we found that the tanh velocity profile with 
small values of 0, which characterize the mean flow close to jet exit, is unstable to a 
large number of azimuthal modes. When 0 increases, the f3 = 1 mode increasingly be- 
comes the dominant instability mode. Further, it is found that the periodic flow is most 
unstable to subharmonic disturbances whose phase velocity is the same as that of the 
fundamental. This matching of the phase velocity only occurs when the jet momentum 
thickness is sufficiently small, since otherwise, the wave speed cannot be matched because 
of dispersiveness. The value 0 = 1/15 is chosen because at this value of 0, not only can 
resonant wave interactions occur, but also the streamwise location implied by this value 
is sufficiently far from the jet nozzle to ensure that a finite-amplitude saturation of the 
fundamental mode can be attained. Some evidence of this particular value of 0 can be 
found in the work of Cohen (1986). 

In Figures (6.8) and (6.9), the subharmonic growth rate and the (relative frequency | for 
various azimuthal numbers are plotted as a function of the amplitude of the fundamental 
mode. The results in these figures are based on Re = 5000, 0 = 1/15, /3\ = 0, T* = 1, 
Mj — 0.1, c *2 = 3.66, = 7.32. These figures show that for a low Mach number jet, 

the overall increment of the subharmonic growth rate of the azimuthal modes j3 = 0, 1 
and 2 increase as A increases. The secondary growth rates of /? = 0 and (3 — 1 are of 
comparable size, while the corresponding quantity for /9 = 2 is relatively small. It is rather 
surprising to note that the threshold amplitude A^ of f3 = 3 mode is rather high, and the 
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phase-locking phenomenon is never established for (3 = 4 mode. The secondary growth 
rate of 0 = 3 is in fact smaller than the growth rate obtained from parallel flow analysis. 
These results are of interest because they indicate that the approximate fulfillment of 
the resonance condition on the phase velocity does not automatically lead to an increase 
in the (temporal) growth rate of a subharmonic mode. Physically, the suppression of 
the growth rate of the subharmonic, and hence the vortex-pairing, means the energy is 
pumped from the subharmonic to the mean flow. It would be interesting to see whether 
these observations could be verified experimentally or by full numerical simulations. 

These results suggest that in order to achieve maximum mixing enhancement in round 
jets via the amalgamation of coherent structures, it is most effective to excite a subhar- 
monic mode with = 0 or ft = 1 modes. 


6.7 Influence of A on the Eigenfunction of Subharmonic 

Modes 

In the previous sections, the secondary growth rate and the |relative frequency! 
in the presence of a base flow consisting of a parallel mean flow and a finite-amplitude 
fundamental mode are discussed. This section considers the influence of the fundamental 
mode on the eigenfunctions of the subharmonic. For clarity, the eigenfunctions of the 
fundamental and its subharmonic are written, respectively, as 

= |fi|expi(^j) (6.4) 

f 2 ,i = |f2,i|expi(0 2 ,i) (6.5) 

f 2 -l = |f 2 ,-l| expt'( 02 ,-i) (6 6) 

where the amplitudes |fi|, |f 2 ,x| and |f 2 _i| as well as the phases 9 U 9 2 ,i and 9 2 _i are 
functions of the radial variable r. 
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For a real & 2 , the substitution of (6.4) - (6.6) into (3.39) and (3.45) yields 

Fx = 2A\f\ | cos(aix + #i) (6.7) 

F 2 = 2e exp(d' 2 ^){|f 2 1 i| cos(a 2 X + /?2<£ + # 2 , 1 ) + |f2,-i| cos(a 2 x - /? 2 <£ + 02,-i)} (6.8) 

Recall that when 02 is purely real, the [relative frequencyj is zero. Further, the eigenfunc- 
tions are normalized in such a manner that A and € becomes the maximum amplitudes of 
the streamwise velocity perturbations. 

In Figures (6.10) - (6.11), the amplitude jf^il and the phase # 2,1 as a function of 
the radial variable, r, are plotted for three different values of A: A = 0.0, A = 0.05 
and A = 0.1. [For reference, the corresponding growth rate and [relative frequency | are 
plotted in Figures (6.2) - (6.3).] As discussed before, these figures characterize the case 
where the resonance condition on phase velocity (6.3) is approximately met, whereas the 
resonance condition on the growth rate is not satisfied. Except for small changes in the 
amplitudes of the radial and azimuthal velocity components as A becomes larger, the 
secondary eigenfunctions are generally unaffected by the size of A. This lends credibility 
to the fact that the mode shape of the secondary instability wave or that of the primary 
during saturation is given reasonably well by the parallel flow results. 

Figures (6.12) - (6.13) show the amplitude and the phase of a near maximally ampli- 
fied subharmonic disturbance (/? = 1) for an unheated low Mach number jet. The jet 
parameters are identical with the previous case. [For reference, the corresponding growth 
rate and “relative frequency” are plotted in Figure (6.4).] These figures represent the 
case where the resonance conditions on phase velocity (6.3) and growth rate (6.4) are 
adequately satisfied. Again, with the exception of a slight decrease of the eigenfunction 
amplitudes in the vicinity of the critical layer as A increases, the amplitude functions of 
the secondary subharmonic modes are minimally affected by the magnitude of A . 

In summary, we conclude that while the finite-amplitude fundamental mode can sub- 
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stantially enhance the growth rate of its subharmonic, it has only a small to moderate 
effect on the mode shape of the subharmonic. 


Phase Velocity , C R Phase Velocity , c, 


Routes of Resonant Interactions: I & II 




Fig. 6.1 Inviscid phase velocity for various mean profiles at 6 — 1/15. 
Solid lines: M, = 0.0, T. = 1. Dot: M, = 0.0, T, = 2. Dash: Mj = 0.8, 
(a) Helical mode: 0=1. (b) Axisymmetric mode: 0 — 0. 
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Comparison of Approaches: Incomprewlble, Isothermal Jet 




Fig. 6.2 (a) Growth rate and (b) |relative frequency! M a function of the amplitude 
of the periodic wave with a near maximally amplified wavenumber. Calculations 
are for = 5000 , 9 = 1/15, Mj = 0.1, T. = 1,/?, = 0, a, = 1.7,/?, = 1. 

Solid lines: multiple scales. Circles: normal mode analysis. 
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Comparison of Approaches: Incompressible, Isothermal Jet 




Fig. 6.3 (a) Growth rate and (b) (relative frequency] as a function of the amplitude 
of the periodic wave with a near maximally amplified wavenumber. Calculations 
are for Re = 5000 , 9 = 1/15, M,- = 0.1, T. = l,/? t = 0, a, = 1 . 7 , 0 2 = 1. 

Solid lines: multiple scales. Circles: normal mode analysis. 





Comparison of Approaches: Incompressible, Hot Jet 
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Amplitude of Periodic Wave , A 


Fig. 6.4 (a) Growth rate and (b) |relative frequency! 35 a function of jthe amplitude 
of the periodic wave with a near neutral wavenumber. Calculations are for 
Re — 500 , 9 = 1/15, M, = 0.1, T, = 2,/?, = 0, a 2 = 3.2, /? 2 = 1. 

Solid lines: multiple scales. Circles: normal mode analysis. 




Compressibility Effect: Isothermal jet 
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Fig. 6.5 (a) Growth rate and (b) |relative frequency! 35 a function of the amplitude 
of the periodic wave with a near neutral wavenumber. Calculations, using normal 
analysis, are for Re = 5000,0 = 1/15, T, = l,/?i = 0. Dot: or 2 = 3.275, /? 2 = 0, My = 
0.8. Dash: a 2 — 3.275, /?j = l,My = 0.8. Chaindot: a 3 = 3.66, /? 2 = 0,My = 0.1. 
Chaindeish: a 2 = 3.66, /? 2 = l,My =0.1. 





Heating Effect: Incompressible Jet 
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Amplitude of Periodic Weve , A 


Fig. 6.6 (a) Growth rate and (b) [relative frequencyl as a function of the amplitude 
of the periodic wave with a near neutral wavenumber. Calculations, using normal 
analysis, are for 0 = 1/15, My = 0.1, 0i = O.Dot: Re = 500, or 2 = 3.2, /? 2 = 0,7’, = 2. 
Dash: Re = 500, a 2 = 3.2, fa — 1 , T, = 2. Chaindot: Re = 5000, a 2 = 3.66, 

/? 2 = 0,T, = 1 . Chaindash: 72e = 5000, a 2 = 3.66, /? 2 = 1,T, = 1 . 



Shear Layer Thickness Effect: Incompressible, Cold Jet 
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Fig. 0.7 (a) Growth rate and (b) (relative frequencyl as a function of the amplitude 
of the periodic wave with a near neutral wavenumber. Calculations, using normal 
analysis, are for Re - 5000, A/y = 0.1, T. = 1,/?, = 0. Dot: a, = 7.3,/?, = 0, 

9 = 1/30. Dash: a, = 7.3,/?, = 1,9 = 1/30. Chaindot: or, = 3.66,/?, = 0,0 = 1/15. 
Chaindash: cr, = 3.66,/?, = 1,9 = 1/15. Solid lines: a, = 1.7,/?, = 1,9 = 1/15. 
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Mode Number Effect: Incompressible, Isothermal Jet 




Fig. 6.8 (a) Growth rate and (b) (relative frequency) as a function of the (unrestr- 
ited) amplitude of the periodic wave with a near neutral wavenumber. Calculations, 
using normal analysis, are for Re = 5000, 6 = 1/15 , Mj = 0.1, T. = l,0 l = 0, 

<*! = 3.66. Solid lines: fa = 0. Dot: fa = 1. Dash: fa = 2. 




Mode Number Effect: Incompressible, Isothermal Jet 




Fig. 6.9 (a) Growth rate and (b) |relative frequency! 33 3 Junction of the (unre- 
stricted) amplitude of the periodic wave with a near neutral wavenumber. Calcu- 
lations, using normal analysis, are for Re = 5000,0 = 1/15, M, = 0.1, T, = 1, 

/?i = 0, c *2 = 3.66. Solid lines: (3 2 = 3. Dot: /3 2 = 4. 
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Effect on |Eigenfunction| 

Incompressible, Isothermal Jet 



Fig. 6.10 Distributions of the (eigenfunction) of a subharmonic disturbance as a 
function of the radial variable r. Calculations are for Re = 5000, T \ = 1.0, My — 0.1, 
0 = 1/15, /?x = 0,a 2 = 1.7,02 = 1. Solid lines: A — 0. Dot: A = 0.05. 

Chaindot: A = 0.1. 




Effect on Phase of Eigenfunction 

Incompressible, Isothermal Jet 
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Fig. 6.11 Distributions of the phase of a subharmonic disturbance as a function 
of the radial variable r. Calculations are for Re = 5000, T. = 1.0, M, - 0.1, 
9 = l/15,/?x = 0,a 2 = 1.7, /? 2 = 1. x: A = 0. •: A = 0.05 . o : A = 0.1. 



Effect on |EigenfuncHon| 

Incompressible, Isothermal .let 
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Fig. 6.12 Distributions of the |eigenfunction| of a near maximally amplified sub- 
harmonic disturbance as a function of the radial variable r. Calculations are for 
Re = 5000, T, = 1.0, A/, = 0.1,0 = 1/15,/?, = 0,a 2 = 3.66, /? 2 = 1. 

Solid lines: A = 0. Dot: A = 0.05. chaindot: A = 0.1. 



Effect on Phase of Eigenfunctions 
Incompressible, Isothermal Jet 
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Fig. 6.13 Distributions of the phase of a near maximally amplified subhar- 
monic disturbance as a function of the radial variable r. Calculations are for 
Re = 5000, T, = 1.0 , Mj = 0.1,0 = 1/15, ^ = 0,ct 2 = 3.66, /3 2 = 1. x: A = 0. 
• : A = 0.05. o : A = 0.1. 
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Chapter 7 

SUMMARY k CONCLUSIONS 


This research concerns the (temporal) instabilities and resonances of the shear layer 
generated by a subsonic heated jet emanating from an axisymmetric nozzle. When the 
shear layer is excited by a harmonically oscillating disturbance of a suitable frequency (or 
wavenumber), an instability wave develops in the flow. For clarity, this instability wave 
is termed a primary wave or fundamental wave. In the initial stages of development — 
where the amplitude of the disturbance relative to the mean flow is small — the stability 
characteristics (e.g., growth rate, phase velocity and mode shape) of the primary wave 
are described quite accurately by linear instability theory. We begin our research with 
the systematic investigation of the linear instability characteristics of a small disturbance 
superimposed on a prescribed parallel mean flow approximated by the hyperbolic tangent 
function. By generalizing Howard’s semi-circle theorem to compressible round jets, we 
have shown that compressibility reduces the range of unstable wavenumbers and has a 
stabilizing effect. The role of subsonic Mach number, My, in affecting jet instability 
properties is numerically examined for the tanh velocity profile with 8 = 1/15. We have 
concentrated our research on this value of 9 because it provides a good representation of 
the velocity profile at a downstream location where the shear layer is sufficiently thin to 
ensure that resonance interactions between two waves (i.e., fundamental and subharmonic) 
can occur. Hence, subsequent discussions, unless otherwise indicated, pertain only to these 
mean profiles with 9 = 1/15. 

Calculations show that the most unstable temporal growth rates of a (2-d) axisymmet- 
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ric mode and a (3-d) first helical mode (i.e., = l)in an incompressible flow are approxi- 

mately equally dominant, and compressibility (in the subsonic range) reduces the growth 
rates of the axisymmetric mode and the helical mode by approximately 19% and 17%, re- 
spectively. Compressibility also decreases slightly the phase velocities of those instability 
modes whose wavenumbers are smaller than the maximally amplified wavenumber, and 
has no apparent effect on the phase velocities of those with higher wavenumbers. In order 
to single out the effect of heating, the stability characteristics of the tanh profile with a 
temperature ratio T«, = 2 — but with Mach number My = 0 — are numerically examined. 
The results indicate that while the phase velocity is rather insensitive to Mach number, it 
is, except for very long axisymmetric waves, substantially decreased by heating. Heating 
also appears to reduce the range of unstable wavenumber and to slightly decrease the 
growth rate of the most amplified wavenumber. The stability behavior of an instability 
wave which has a large wavelength relative to the radius of the jet nozzle is compared to 
that obtained from a heated vortex sheet analysis. The results indicate that the phase 
velocities of long waves (i.e., a — ► 0) — irrespective of the azimuthal mode number — 
in a tanh velocity profile agreed well with those obtained from the vortex sheet. For the 
corresponding growth rates, this agreement is found only in axisymmetric modes but not 
in helical modes. In fact, the disagreement becomes progressively worse as the helical 
mode number and the shear layer thickness is increased. 

When the primary wave propagates downstream, its amplitude will — as a result of 
the exponential growth rate predicted by the linear instability theory — quickly grow to 
an appreciable size. Because of the mean flow divergence (caused by viscous spreading) 
and nonlinear effects, this wave eventually reaches a finite-amplitude equilibrium state 
(i.e., saturation) at some downstream location. We approximate the nearly periodic flow 
which arises from the saturation of the primary disturbance by a parallel mean flow of 
the hyperbolic-tangent type and a periodic component whose phase velocity, wavenumber, 
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mode number and mode shape are the same as those associated with the primary wave. 
The principal task of this research is to determine the resonances and instabilities of this 
periodic flow. The conditions for resonant interaction between two arbitrary instability 
modes are derived. The resonant interaction between two modes allows an effective mech- 
anism for the transfer of energy from the mean flow to an instability mode. In order 
to satisfy the resonance conditions, we have shown that the wavenumbers and the mode 
numbers of the two unstable modes, a primary and a secondary, must be restricted to the 
following choices: 

1. The primary mode is axisymmetric and has an axial wavenumber close to a neutral 
wavenumber. 

2. The secondary wave, which can have an arbitrary mode number, is almost the max- 
imally amplified subharmonic wave; the wavelength of the subharmonic is twice that 
of the primary wave. 

A key result from the resonant conditions indicates that resonant interaction between two 
modes can occur only when the jet momentum thickness is sufficiently small since, other- 
wise, the wave speeds cannot be matched because of dispersiveness. When the resonance 
conditions are met, a subharmonic mode interacts strongly with the primary to produce 
an instability mode whose wavenumber and mode number is the sum of the wavenumbers 
and mode numbers of the two interacting waves. If the subharmonic mode is axisymmet- 
ric, it will interact with the fundamental to reproduce itself. The interaction involving a 
helical mode with mode number f3 will, however, produce a mode with mode number 
and vice versa. Therefore, whenever an interaction involves a subharmonic mode which is 
helical, we shall assume that the secondary disturbance is comprised of a pair of helical 
waves whose propagating angles are equal and opposite. As a result, the secondary distur- 
bance, whether axisymmetric or helical, when interacting with an axisymmetric primary 
wave will reproduce itself and therefore increase its growth rate. We have obtained this 
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(modified) growth rate of the subharmonic as a function of the amplitude of the primary 
wave — for different jet Mach numbers and jet temperatures — using two independent 
methods; the method of multiple scales and a generalization of the normal mode analysis 
in a spatially periodic base flow. 

We found that when the resonance conditions are met almost exactly and when the 
amplitude A of the primary mode is less than 3 % of the mean flow, the two approaches 
are in good agreement. However, as A increases from 3 % to 10%, the perturbation 
method appears to progressively over-estimate the subharmonic growth rate. It is also 
shown that when the amplitude of the primary increases from 0 to 2%, the growth rate 
of the subharmonic mode is augmented by about 7%. This leads us to believe that the 
secondary instability — as manifested by the increase of the growth rate of the subhar- 
monic — will enhance the pairing of the vortex rings associated with the primary wave. 
Our results indicate that the secondary (subharmonic) growth rates for a (3-d) instability 
wave with azimuthal number f3 = 1 and a (2-d) axisymmetric wave are equally dominant, 
while those associated with higher spinning numbers are relatively insignificant. Further- 
more, compressibility has a moderate stabilizing effect on the secondary instability modes 
(whether axisymmetric or helical) so that vortex pairing, which results from subharmonic 
instability, is expected to be less pronounced in compressible flows. 

In order to have a representative idea of the effects of compressibility and heating, 
we use the total increment, denoted by N, of the growth rate of a subharmonic mode 
when the amplitude of the primary wave changes from 0 to 10%. Compressibility reduces 
the increment, N, of an axisymmetric mode and a helical (/? = 1) mode by about 15% 
and 12%, respectively. While heating seems to suppress moderately the increment, K, 
of the secondary growth rate of an axisymmetric mode (by about 15%), it reduces the 
corresponding increment of the /3 — 1 mode significantly (by about 40%). This leads 
us to conjecture that heating tends to inhibit vortex-pairing in an axisymmetric jet, and 
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this inhibition is most pronounced for the /? = 1 mode; we infer that vortex-paring in a 
heated, incompressible jet is most likely to be a (2-d) axisymmetric phenomenon. Our 
calculations also show that the mode shapes of the subharmonic instability modes are 
generally unaffected by the presence of the (finite) amplitude primary wave. Finally, 
when we calculated secondary instability for other values of 0, we found that the growth 
rate of the subharmonic mode increased in proportion to 1/0. This implies that vortex- 
paring is more likely to occur near the jet nozzle, (of course, we also have the requirement 
of phase speed matching which is most perfect for small 0). 

In summary, the principal effect of excitation is to generate a saturated large-scale 
structure whose vortices may merge because of subharmonic (secondary) instabilities. 
Since the secondary growth rates of the subharmonic instability modes — which are at 
least partially responsible for jet spreading and mixing — are suppressed by heating and 
compressibility, we conjecture that it will be more difficult to control vortex-mergings in 
a subsonic heated jet than in a low speed cold jet. 
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APPENDIX A 


The Primary Stability Equations 


The stability equations for the evolution of linear disturbances superimposed on a 
parallel mean flow are of the form 

Ajfj = 0 (2.31) 

where Aj is a 5 x 5 matrix whose elements, a,y, are 
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the primes denote differentiation with respect to r, and the subscript 0 refers to the parallel 
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APPENDIX B 

Stability Equations for the Primary Eigenmodes 

The elements of Ai and Bj in equation (2.33) are 
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APPENDIX C 

Governing Equations for the Secondary Disturbance 

Perturbation Method 

The governing equation for the secondary (subharmonic) disturbance — obtained at O (5e) 
— is of the form 

I»lj(P)l 2 = J5+R+ + E-H.- + c.c. (3-4) 

where 

R+ = [F u , F v , Fxu, Fp, Fp\ , 

R— = [F u — , F v — , F w —, Fp- , Fp~] , 

and the components of R+ and R_ are found to be: 
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In the above equations, the primes denote differentiation with respect to r, and the 
tilde represents the complex conjugate of a quantity. 
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APPENDIX D 


Amplitude Equation Derived from Viscous Analysis 


Here we consider the solvability condition for a system of ODE’s. The approach here 
follows closely from Nayfeh ( 1980 ). In analyzing the stability problem of a compressible 
and viscous jets, the lowest inhomogeneous problem occurs at order This problem 
governs the quadratic interaction between two linear waves: Fi and F2. The appropriate 
governing equations for the case where the interaction between Fi and F2 reinforces F2 
are 
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where F U -,F V -,F W - and S p _ are given in Appendix E. 

The boundary conditions at r = 0 are 

Ui 2 (0) = Pl2(0) = t>'i 2 (0) = “'nC 0 ) f° r & = (^- 5 ) 

u'i 2 (0) - p'i2(°) = U 12(0) = W12(°) for odd (S.6) 

and at r — ► oo, the disturbances vanish. The inhomogeneous problem has a solution 
provided that a solvability condition is satisfied. To determine this condition, we multiply 
S.l by r£i, S.2 by rf 2 , D. 3 by rf 3 and DA by r£ 4 . After adding the resulting equations 
and integrating the results by parts, we obtain 
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The adjoint equations are defined by setting each of the coefficients of ui 2 , iq 2 , and 
pi 2 in the integrands of D. 7 equals to zero. They are 
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+ ~ C2Po ~ fler*)^ 3 ~ - ^' rp °^ 4 = 0 (Z?.10) 
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To determine the boundary conditions on the £'s, we set the boundary terms in D,7 = 
0, and after invoking the homogeneous boundary conditions (D.5) and (D.6), the adjoint 
boundary conditions are found to be 


6(0) = 6(0) - 6(0) = 6(0) = 0 for ($2 = odd (D. 12) 

6(0) = 6(0) = 6(0) = 6(0) = 0 for ft # odd (D.13) 


It may be noted that if we replace 6 by t*i 2 , 6 by -V12, 6 by W12 and 7P06 by pn in 
(D.8) - (D.12), the asymptotic equations where contant flows prevail are identical with 
(D.l) - (D.4). Comparing the adjoint problem with the original homogeneous problem, 
we find that although their boundary conditions are identical, the differential equations 
are different, and hence, the homogeneous problem is not self-adjoint. The solvability 
condition, which obtained from equation (D.7), is given as 
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where the superscript v indicates that a quantity is obtained from the viscous analysis, and 
the tilde denotes the complex conjugate. Here gj, gj an< ^ & Tlse from the dependence 
of the F 2 on the slow scales, gj contains the particular set of interactions between Fi 
and F 2 that reproduces F 2 , while gj, g| are included in the above equation to allow the 
resonance to be developed spatially; their explicit expressions are not given here because 
our analysis is based on temporal theory. 

Again, as in equation (3.32), the same exponential dependence of B v on the slow scales 
is assumed, i.e., 


B v ~ exp(Ajxi + Ajti + Aj^i) 


The above equation, together with (ZX15) and its complex conjugate, yields 


(£U6) 




’AS’+kll 


l Af+(2j;5?)sAliAS+(2««S)i.A;AS+(25; 9 S s A;A5 = |,;| ! (£>.17) 


The A u ’s are in general complex and ( • )# denotes the real part of a complex number. 
We consider a subharmonic disturbance with fixed wavenumber and it is allowed to grow 
only in time, (D.17) can then be simplified to 


A 


V 

2 



{D. 18) 


The physical interpretation of A 2 has been discussed in Chapter 3, and will not be repeated 
here. 



174 


APPENDIX E 


Resonant Wave Interactions Terms 


The governing equation for the linearized disturbance equation for pressure — 
obtained at O(Se) — is of the form 


t ( \ _ ^Pl2 

MP12) — ~ ~ — + 


dr 2 


r «o ~ Cl 


P 0 
Pol 


dp 12 
dr 


+ 


where 


a^(uo - C 2 ) 2 p 0 M* ~ a] 


dB dB dB 

rh 8 = + ^2^ h ^stt - + ^12-B 

atx v2?i C/0i 


P12 = r/is 

( 3 . 29 ) 

( 3 . 30 ) 


A 1 214* (j[ 

hn = ta 2 F u - + [ — + — + i/3 2 F w -'/r — ia 2 (uo — C 2 )poM 2 F p - 

r uq — O 2 po dr 

In the above equations, all variables, except for F u _, F v -, F w - and F p _, have been 
defined in Chapter 3. The undefined variables Me 

F u - = -{po[v-2«x + viu'_ 2 + 

+t/9 2 U-2 7*“ + «OilU_2Ul - IQ!2«l«-2] + «o(w-2Pl + ViP- 2 ) 
+U 0 (P-2*O'l«l ~ Pl«0!2«-2) + Pl<^2«-2 + P-2«l^l) 

F v - = — {/>o[v— 2^1 + + tCKiU_2Vl 

— ia!2 u l v -2 + tPi—r 1 + *-%rV- 2 wi — 2u '~ 2U ' 1 . ] 

+ti 0 (P- 2 ‘Q!iVi - piia 2 v- 2 ) + p\o 2 v- 2 + p- 2 via\} 
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Fw- = -{po[v~2V)\ + v\w'_2 + «aitt-2U>i — »a2tiitu_2 

+ifil + ifi 2 2^ + =^2l + Si2=l] 

+Uo(p-2*C*lU>i " P\i&2W —i) + P-2 <t 1 w 1 + P l a 2^-2{ 

Fp _ = -{v_2pi + Vip'_2 + u_2»aipi - Uip_2»a2 
+ ^P1»^1 + 3 7 L p-2*^2 + 7P-2(^ + u i+ 

+ Uitai) + 7 Pi(^t 2 - + v!_2 + $*>-2*02 _ «-2* a 2)} 

In accordance with the discussion in Chapter 3.2, the eigenfunctions u-2> f-2> 2»P— 2 

and p - 2 are for an unstable mode with a wavenumber —02 and a mode number @2 ■ 
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APPENDIX F 


Secondary Instability Equations Using Normal Mode 

Analysis 


The shape function of the subharmonic disturbance g is written as 

g = ( U 2 ,U-2,V2,V-2,W2,W- 2 ,P2,P-2,P2,P-2) T (3.47) 

The stability of subharmonic disturbance superimposed on a periodic flow consisting of a 
finite amplitude primary wave and a parallel mean flow is governed by the following set 
of 10 coupled linear ODE’s: 


§U2 — Pqu! q V 2 - tOL2P2 — AFu- = 5^2(poU2 + Apitl- 2 ) 


(F.l) 


1 2»ft 

(§ - - ']fa2 W Z - D P 2 - ~ ^2(P0«2 + ApiU_ 2 ) 


{F- 3 ) 


+ (§ - & 2 P 0 - ^~2') U, 2 - ~~~P2 - AFyi- = £ 2 (po«>2 + ApiW-i) (F.5) 


-7P0«a2«2 - 7Po(^ + D) V2 - 7P °^ 2 W2 + a 4 4p 2 + a<sP2 - A(F P _ + V’con) = ^ 2 (p2 + Api/po) 


{F.l) 


I / 1 , Poi02 . . f, 

-po ,Q! 2 u 2 - ^o y 2 - Po\ — + D) v 1 — u>2 - «0*<*2/>2 - AFp- = <J 2 P2 


§ = ^ - a 2^j - P0*’ot2«0, 


(F. 9) 


where 
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After replacing the subscript 2 of the components of the shape function g by -2 (and 
-2 by 2) and by substituting the coefficients (with the exception of tfa) by their complex 
conjugates in equations (F.l), (F.3), (F.5), (F.7), and (F.9), equations (F.2), (F.4), (F.6), 
(F.8), and (F.10) are obtained, respectively. 

In equations (F.l) - (F.10), the prime denotes a differentiation with respect to r, uo 
refers to the mean velocity in a moving frame of reference, <144 and <145 are given in 
Appendix B, P u ->F v ->F w - } F p - and are given in Appendix E, and tp cor ly which arises 
from the quadratic interactions in the conduction term of the energy equation, is found 
to be 


'Peon = -^Z'Pj 

3 = 1 


(F- 11) 


where 


, Piot\p-2 P\ot\p~2 4 p 0 Q!ia 2 pip _ 2 , 2 a 1 a 2 pip _ 2 2 , P0P1P-2 , P1P-2 

= — 7f 7o 7 + 7 + 

PoPlP-2 4PlPoP'-2 , Pi P-7 , 2 p' 1 p '_ 2 PQp\PlP-1 4 poPjp !_ 2 
Va = -TZjSx 3 + 773T + 


Vp I) 

V>3 = - 


P 0 


M) 


Po 


(rVo) 


Po 


Pl<#2P-2 4poPlP-2^2^1 2/?i/? 2 p_ 2 pi _ PoQiiPlP-2 


(r 2 P 2 o) 


(r 2 P 3 0 ) 


(r 2 P 2 o) 


Po 


'Pi = 


_ POPiP-2 1 PoPjP-2 PO/9 2 PlP-2 PlPpP-2 PlPpP-2 . 2PxPqP-2 tyjP-i . P-2 


pT 


(••Po) ( r2 Po) 


pi 


( r Po) 


(Po) W ( f2 Po) 


V’S = 


2a\pip~2 2p"p_ 2 P- 2 <*iPi , 2p_ 2 pia l a 2 f P- 2 P 1 4p_ 2 poPx _ P- 2 P 1 , 2 P > _ 2 Pi 


Pi 


p! 


1 — + 

Po 


?o 


+ — r 
Po 


pT 


M) 


Po 


_ p-2piPi | 2^ 2 /3iP-2Pi , P- 2P0P1 t P-2 PqPi , 2P2P0P1 2p_ 2 pi | 2/3jp_ 2 pi 


V'6 = - i 


( f2 Po) 


( r2 Po) 


+ 


Po 


+ 


( r Po) 


+ 


Po 


( r Po) + ( r2 Po) 


, 2o! 2 p_ 2 pi 2p"jPi 

* = — ^ ?r 


The appropriate boundary conditions at r — ► 0 are 

U 2 (0) = p 2 (0) = p 2 (0) = V 2 (0) = «>2( 0 ) for & = odd 
ti ' 2 (0) = p' 2 (0) = p' 2 (0) = v 2 (0) = to 2 (0) for P odd 


(F. 12) 
(F.13) 
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u_ 2 (0) = P— 2(0) = P— 2 (O) = t/l 2 (0) — u;L 2 (0) f° r Pi ” °dd (F.14) 

u'_ 2 (0) = P- 2 (0) = p 1 - 2 (0) = t>_ 2 (0) = u;_ 2 (0) for ft 7^ odd (F.15) 

and at r — ► 00, the subharmonic disturbances vanish. 
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